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Synopsis 


From the very early days of high-T c superconductivity it has been believed that sub- 
stitutional studies in cuprates can help in understanding superconducting and magnetic 
properties of these materials [1]. In view of this, a number of experimental studies have been 
done on these systems doped with impurities - both magnetic, e.g. Fe +3 , Ni +2 , Co +3 and 
nonmagnetic, e.g. Zn +2 , Al +3 , Ga +3 - that substitute copper in the copper oxide planes. 
Some of the most interesting effects of doping these materials is that irrespective of the mag- 
netic nature of the dopant, local moments are formed in the system, T c is drastically reduced, 
and nonmagnetic impurities reduce Neel temperature TV. As the parent compounds of all 
these materials are antiferromagnets (AF) and the antiferromagnetic correlations persist 
even in the superconducting phase, formation of local moments in presence of nonmagnetic 
impurities can be naively understood as due to removal of copper spins from a compensated 
spin system. But this picture of AF plus spin-hole becomes simplistic in view of the recent 
NMR measurements by Mahajan et al. [2] which suggest that the moment resides mainly 
on the four copper sites neighbouring the impurity. So a more detailed picture is called for. 
In this thesis we try to develop a deeper understanding of the formation of moment and to 
gain insight into other effects of impurity-doping in high-T c cuprates. 

The thesis begins with Chapter 1 giving a general introduction to the high-T c cuprates' 
and to the particular questions about the impurity-doped systems that we address. This 
leads us to the problem of nonmagnetic impurities in a Mott-Hubbard AF in Chapter 2. A 
formalism for Hartree-Fock (HF) mean-field treatment of the Hubbard AF had already been 



loped [3]. We study the effects of impurities in this state within a T-matrix analysis. It is 
d that such impurities give rise to energy-states within the Hubbard gap (defect states), 
ng amplitude mainly on the four sites neighbouring the impurity. These, if occupied by 
irons, give rise to local moments as seen in experiments. Moments are seen to form even 
ie superconducting and ‘metallic’ phases when there is no AF long-range-order (LRO) 
re system (although short-range-order persists). So in Chapter 3 we address the problem 
nonmagnetic impurity in a small cluster of Hubbard model having short-range-order, 
find that the existence of the defect state is robust with respect to the range of order in 
system. 

The next question to ask is how these impurities affect the magnetic dynamics in the 
system and thereby the Neel temperature. Chapter 4 of the thesis deals with impurity- 
:tering of spin waves in an impurity-doped Hubbard AF. In a two-dimensional (2D) system 
find a logarithmically diverging correction to the long wavelength magnon energies. In 
“ of an anisotropic (layered) 3D system this divergence is cutoff by the interlayer hopping, 
use linear spin-wave theory to calculate the reduction in sublattice magnetisation due 
:hermal excitation of spin waves at finite temperature and find the dependence of T N on 
>urity concentration, which agrees well with neutron-scattering experiments of Keimer et 
[4]- 

If a hole is added to the system having an impurity, it will go in the defect state as that is 
: lowest ener gy state available for a hole and the moment will be lost. But fluctuations can 
nsfer some of the spectral weight of the hole from the defect state to the upper Hubbard 
id so that part of the moment can be recovered. In Chapter 5 we study how the local 
ment is recovered due to thermal and quantum fluctuations. 

Few questions still remain unanswered. For example, how the moment survives at low 
aperatures in presence of impurities, why do all nonmagnetic impurities, irrespective of 
?ir electromc structure, produce the same moment and why is this moment less than 
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what one would expect for an 5 = ^ spin-hole in an S = ~ AF. We try to answer these by 
studying the effects of nonmagnetic impurities in a charge-transfer insulator(CTI). Chapter 6 
of the thesis contains the results of our self-consistent numerical HF calculations on pure 
and impurity-doped CTI described by a three-band Hubbard model. In absence of direct 
oxygen-oxygen hopping (tpp), the system shows nesting property and goes into an AF state 
for any non-zero positive Coulomb repulsion on the copper sites (U d ). The nesting, however, 
is lost in presence of a finite t pp and the system becomes AF only above a certain critical 
value of Ud, which increases with t pp . The transverse spin excitations in the system have a 
gapless Goldstone mode. Each nonmagnetic impurity in a CTI produces two defect states. 
The upper defect state remains occupied by electrons even in presence of finite concentration 
of mobile holes and the local moment is robust. Density of the occupied defect state on the 
copper sites is less than 1 and this density decreases for a larger copper-oxygen hopping. 
This clearly shows that part of the local moment escapes to the oxygen sites because of 
Cu3d-02p hybridisation. The local magnitude of the local moment is also found to be quite 
insensitive to the impurity potential V suggesting that the observed moment is independent 
of the detailed nature of the impurity. 

In Chapter 7 we develop a formalism to treat magnetic impurities in a Mott-Hubbard 
AF describing 3d electrons in the copper-oxide planes of high-T c superconductors. In the 
case the impurity spin (S') equal to the host spin ( S ), represented via modification in the 
local hopping strength around the impurity, we get a momentum independent non-diverging 
renormalisation of the ntagnon energies. In the long wavelength limit, our result agrees 
with that of Wan et al. [5] on a Heisenberg AF. We also present our study - using exact 
eigenstates method - on the effects of magnetic impurities on magnon wave functions and 
density of states. When S' # 5, we develop a representation of such impurities in terms of 
multiple orbitals on the impurity site and examine their effects on the magnon spectrum via 
a scaling study. 
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Chapter 1 


Introduction 


The revolutionising discovery of high-temperature superconductivity in cuprates by Bed- 
norz and Muller [6] generated tremendous amount of interest and enthusiasm, and attracted 
the attention of physicists and material scientists from all over the world. So powerful was 
the wave that the first few years of high-T c superconductivity saw tens of thousands of 
publications on the subject. But the period of initial euphoria is over and the task of under- 
standing these materials is proving to be more and more difficult and challenging. Ten years 
after Bednorz and Muller’s discovery, it seems that although some of the basic properties of 
these compounds have been unambiguously ascertained, there is till some time before the 
mechanism of superconductivity is understood properly. 

I 

One of the many techniques used to study high-T c cuprates is substitution of Cu in the 

I 

I 

Cu0 2 planes of these materials by different impurities like Fe, Ni, Co, Zn, Al, Ga etc. Even j 
nonmagnetic impurities like Zn +2 , Al +3 , Ga +3 induce local moments in these systems. They 

j 

also have drastic effects on superconducting transition temperature (T c ) and reduce Neel 1 
temperature (Tn) of these materials. Some theoretical studies have been done to understand 
the magnetic effects of these impurities. For example, Wan et al. have done a detailed 
study on the effects of impurities on the magnetic dynamics in a Heisenberg antiferromagnet 
(AF). But to our knowledge, no attempt has been made to understand the formation of 
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il moments because of nonmagnetic impurities in the Cu0 2 planes, how they affect Tn 
why all nonmagnetic impurities, irrespective of their electronic structure, give the same 
ment. In this thesis we address these key questions and other related problems involving 
)urity-doping in high-T c cuprates within an electronic model such as a one-band or a 
ee-band Hubbard model wherein electronic and magnetic properties can be examined on 
equal footing. 

Before w r e go on to discuss our work, we give a brief introduction to the class of materials 
are interested in, some of their basic properties and the particular questions we are 
dressing. 

.1 Structure 

The first high-T c cuprate compound discovered by Bednorz and Muller, La 2 _ x Ba x Cu04, 
s the quasi two dimensional (2D) structure of K2N1F4 (Figure 1.1). It has a series of Cu0 2 
anes with charge reservoirs in between (LaO layers in this case). The parent compound 
i 2 Cu0 4 and also other derivatives, such as La 2 _ x Sr x Cu0 4 , possess this structure with elon- 
tted Cu-0 octahedra. The Cu-0 distance is rather short, ~ 1.9A. La 2 Cu0 4 is orthorhombic 
room temperature and becomes tetragonal at ~ 500K. The temperature at which this 
ructural phase transition from an orthorhombic to a tetragonal phase occurs depends on 
ie doping concentration x in La 2 _ x M x Cu0 4 (M=Sr, Ba) [7]. 

The superconducting Y-Ba-Cu-0 also has a series of Cu0 2 sheets in the ‘ab’ plane (as 
is usually referred to in literature). There are two sheets per unit cell, ~ 3.2Aapart. These 
vo planes are 8.2Aapart from the next two planes. In addition it has Cu-0 chains along 
ie b-axis. The orthorhombic lattice parameters of YBa 2 Cu 3 0 6+ ,s (1-2-3) vary with oxygen 
mtent and when 6 ~ 0.6A, the structure becomes tetragonal. In YBa 2 Cu 3 06 (S = 0.0) 
iere are no chain oxygens and Cu in the chains are in the +1 state. Structural studies 
uggest that oxygen vacancies may be ordered when 5 = 0.5 and 8 = 0.25. In 6 = 0.5 
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, where A=Ba, Sr, Ca, ... Figure taken from 
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position, fully oxygenated CuO chains are present along the b-axis alternately with fully 
Tgn deficient chains. Because of orthorhombic structure these compounds show extensive 
ming. Across the twin boundary there is a 90° rotation of the a- and b-axes. 

All the Bi- and Tl- based high-T c compounds discovered later possess Cu0 2 planes as an 
rntial structural unit. These Cu0 2 planes dominate the high-T c material from magnetic 
l electrical points of view and also, to a large extent, from the structural point of view. It 
;enerally accepted that superconductivity in these materials occur in the Cu0 2 planes and 
ood understanding of the Cu0 2 planes is required for an explanation of superconductivity 
these materials. Since the in-plane and out-of-plane structure, Cu-0 distance etc. are 
istically different, one can expect these materials to be anisotropic. In fact, the ab plane 
d c-axis (perpendicular to the ab plane) properties in the compound are, in some cases, 
alitatively different (see discussions below). 

.2 Phase diagram 

The parent compounds of the cuprates are ‘doped’ which ultimately leads to metallic 
haviour and superconductivity. Doping is achieved either by heterovalent substitution as 
La 2 _ x Sr x Cu04, in Nd 2 - x Ce x Cu0 4 and in Bi 2 Sr 2 Ca 1 _ x Y x Cu 2 0g or by a variation of the 
tal oxygen content as in YBa 2 Cu306+,5. Doping introduces additional charge carriers into 
e Cu0 2 planes [8]. In La 2 _ x Sr x Cu04 the formal valency of strontium is +2. To maintain 
large neutrality in this case, the in-plane oxygen atoms change from an 02p 6 to an 02p 5 
ate leaving additional holes in the planes. In Nd 2 _ x Ce x Cu0 4 the neodymium Nd +3 state 
replaced by a cerium Ce +4 yielding an electron which enters the planes. The sign of 
le measured Hall coefiicient supports the picture of hole carriers in La 2 _ x Sr x Cu0 4 and of 
ectron carriers in Nd 2 _ x CexCu0 4 [9]. 

In YBa 2 Cu306+i doping amounts to the addition of oxygen atoms in between the out-of- 
lane Cu atoms. This leads to the formation of Cu-0 chains. It is believed that these oxygen 
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Metallic «— ■ ■ Insulating Metallic 



Figure 1.2: Temperature vs. doping concentration phase diagram for YBa2Cu 3 0 7 _5 (left 
panel) and of La 2 _xSr x Cu0 4 combined with Nd 2 _ x Ce x Cu0 4 (right panel). Figure taken 
from Reference 46. 

atoms form 0~ 2 states by absorbing electron from the rest of the material and, to a certain 
extent, from out of the Cu0 2 planes. This process is controversial. If true, it introduces 
the possibility of hole type carriers in the Cu0 2 planes and it is certainly responsible for 
the non-linear relation between the dopant concentration and the in-plane carrier density in 
■YBa 2 Cu3C>6+<5 [10]. 

In Figure 1.2 we have shown the phase diagrams of La 2 _ x Sr x Cu0 4 , Nd 2 _ x Ce x Cu0 4 and 
of YBa 2 Cu 3 0 6+< j in the temperature (T) vs. doping concentration (x) plane. Depending on T 
and x the cuprate superconductors can be varied continuously between an antiferromagnetic 
insulator, a spin-glass phase, a high-temperature superconductor and a ‘normal’ metal. In 
La 2 _ x Sr x Cu0 4 , for example, upon doping, antiferromagnetism vanishes rather rapidly (~ 2% 
of hole). It becomes superconducting at about ~ 6% of hole doping. In the intermediate- 
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ing range, it shows a spin-glass behaviour. Above T c it behaves like a ‘strongly correlated 
,al’. Upon further doping more standard metallic behaviour shows up [10]. 

3 Properties of high-T c cuprates 

3.1 Superconducting properties 

These copper-oxide superconductors share all the basic properties of ordinary super- 
lductors like Meissner effect, zero resistance and pairing of quasiparticles with charge 2e. 
iring is seen in Shapiro steps displayed by Josephson tunnel junctions biased with both 
and dc voltages. These steps have voltage spacing hv/ 2e which is expected if quasiparti- 
s of charge e are bound into pairs [11, 12, 13]. Also, the flux-quantum has been directly 
iasured and is hc/2e, indicating pairing [14, 15]. Moreover -Andreev scattering experiment 
zeals pairing of particles with opposite momenta. Also, Josephson tunneling has been 
served to occur between a Pb-Sn alloy and YBa2Cu 3 0 6 +5 [16]. Such an effect is thought 
be impossible between a singlet-paired and a triplet-paired superconductor [17]. Since, a 
)-Sn alloy is singlet-paired, YBa 2 Cu 30 6 +,5 is also thought to be singlet-paired. 

In measurements of electron tunneling and electromagnetic absorption, an energy gap 
s been observed in the distribution of the energy levels available to the system, and it 
of the same order of magnitude at temperatures T <CT C as that predicted by the BCS 
eory. The observed values initially extended over a range from about 2 to 14 k B T c , later 
periments put the value in the range 4 - 8k B T c . Recently there have come mounting 
idence for an anisotropic gap and a d-wave pairing in the high-T c superconductors [18]. 

.3.2 Some unusual properties 

The high-T c superconductors have three most unusual fundamental properties — large 
: , short coherence lengths and large spatial anisotropy. The large T c values open the 
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way for many new applications which were not previously attractive economically. From a 
fundamental viewpoint, T c values of the order of 90K or higher are probably too high to 
be explained by BCS theory unless one assumes that some interaction much stronger that 
electron-phonon interaction is at work [19]. When T c is large, many types of excitations, 
such as phonons, are present in the upper temperature ranges of superconductivity. The 
presence of these excitations, if they break Cooper-pairs, is expected to affect some of the 
properties of these compounds, such as, T c and the critical current density J c [20, 21, 22]. 

In these new materials £gl, the Ginzberg-Landau coherence length at T <C T c is typically 
between 0.5 and 30A, depending on the crystallographic direction of the electron momentum, 
the substance, the type of experiment etc. [7]. In the conventional superconductors, £gl( 0) 
is typically a few thousand Angstroms, unless it is decreased by alloying or defects. Because 
€gl(T) is much shorter than electromagnetic penetration depth A in high-T c materials, they 
are all type II superconductors with extremely high values of the upper critical field H c2 . 
Also associated with the short coherence length is a weak pinning of fluxoids in the mixed 
phase . [23], compared with pinning in the conventional superconductors. This decreased 
flux-pinning diminishes the size of J c . 

The high-T c cuprates have a high spatial anisotropy. For Y-Ba-Cu-0 the least anisotropic 
among the high-T c superconductors, one finds = 14±2A, <f c = 1.4±3A, X ab — 1400Aand 
A c ~ 7000 A. In the ‘normal state’, the electrical resistivity data show similar anisotropy. In 
Y-Ba-Cu-O, resistivity in the c-direction, p c , is typically larger by a factor of 50 or more 
than the planar resistivity p Qi) . For a Bi-Sr-Ca-Cu-0 sample (with T c ~ 81K) this anisotropy 
is extremely large, nearly 10 5 [24]. For homogeneous single-phase samples which are super- 
conducting, the following general differences are seen — a) p ab is ‘metallic’ with ^ positive 
and comparable to b) p c is usually ‘non-metallic’ < 0) [24], 

One of the basic questions about the normal state of the high-T c cuprates is whether 
they are standard Fermi liquids or not. On the one hand, there is a well-defined Fermi 
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rface from photo-emission [25] and positron-annihilation [26] experiments, and the effective 
iss determined from the magnetic susceptibility [27] and the Drude component of optical 
nductivity [28] are consistent with each other. On the one hand, the transport properties, 
:e resistivity and Hall effect, can be interpreted only in terms of a doped Mott-insulator, 
ther than a standard Fermi liquid. How to reconcile these two apparently contradictory 
pects is the main challenge to theory. 

The temperature dependence of resistivity in several families of high-T c superconductors 
ow the following general features it increases linearly over a wide range of temperature, 
is quite small when extrapolated to absolute zero of temperature. In conventional metals 
e temperature dependence of resistivity can be divided into two distinct regimes. In the 
gh temperature regime, calculation of resistivity shows p — A + BT for T > 0.2 Q D , where 
d is the Debye temperature. At low temperatures one gets, pT 5 , for T <C 0.2Q D . Extension 
the linear temperature dependence of p for Bi-Sr-Ca-Cu-0 down to 10K suggests that some 
echanism other than electron-phonon collision is at work. 

The behaviour of the Hall constant is even more peculiar. Figure 1.3 shows R H changes 
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with composition in electron and hole doped superconductors. The signs of Hall coefficient 
indicates that in the La-system the carriers are holes and in Nd-system they are electrons. 
This observation is consistent with a doped Mott-Hubbard picture. The magnitude of Rg 
has a i dependence, but deviates from this at higher doping concentrations [29]. The most 
interesting point is the change in Rg-sign at a critical concentration of doping where super- 
conductivity disappears and usual metallic behaviour shows up. The region of ‘anomalous’ 
Hall constant coincides with the regime of linear temperature resistivity, as well as with 
linear temperature dependence of R^ 1 [30]. 

The angle-resolved photoemission experiment is a useful tool to study the energy spec- 
trum of electrons. From the energy and angular dependence of emitted electrons one can, in 
principle, recover the spectral function A(k,uj). First in Bi-Sr-Ca-Cu-0 [25], then in Y-Ba- 
Cu-0 [31] compounds a well-defined Fermi surface has been found in good agreement with 
local density functional calculations. This Fermi surface is large as required by the Luttinger 
theorem. The more subtle question is the precise shape of the spectral function which can 
check the Fermi liquid theory in a quantitative way. Unfortunately, the present precision 
does not allow one to make unambiguous conclusions. The main puzzle comes from the 
comparison with transport measurements. As mentioned above, they are consistent with a 
Mott-Hubbard picture, ie., the charge carrier concentration is x instead of 1—x, where x is 
the doping concentration. In other words, the Fermi surface, even if it exists, should be a 
‘small’ one corresponding to these ‘pockets’. 

1.3.3 Magnetic properties 

From neutron scattering studies on La 2 Cu0 4 it was concluded that the system under- 
goes a three- dimens ional antiferromagnetic transition at the Neel temperature T N , at about 
200K. The staggered order parameter has a characteristic temperature dependence of a 3D 
transition. This transition was first discovered by Vaknin et al. [32] in their powder diffrac- 
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Lq 2 - x Sr x Cu 04 



gure 1.4: Correlation length vs. doping concentration for La2- x Sr x Cu0 4 . The solid curve 
calculation of the mean spacing between the holes. Figure taken from Reference 38. 

>n study. It was found that the 3D transition temperature is extremely sensitive to oxygen 
ncentration [33]. Although the maximum possible Tn may be ~ 300K, a small change in 
:ygen concentration can reduce Tn drastically. Vaknin et al. [32] also found that in the 
nit T — > 0 the ordered moment is 0.5 ± 0.15 /j-b, considerably smaller than the saturation 
iue for Cu +2 , which is ~ 1 p#. 

Y-Ba-Cu-0 systems show striking similarities with lanthanum La2_ x Sr x Cu0 4 system, 
be first evidence of magnetic ordering in Y-Ba-Cu-0 came from muon spin rotation (/iSR) 
periments by Nishida et al. [34]. Soon after that superlattice reflections indicative of 
D long-range antiferromagnetic order (AFLRO) were observed by Tranquada et al. [35] in 
iutron powder diffraction measurements. The spin structure which was deduced from the 
:ak intensities was later corroborated by several groups [36, 37], 

Addition of holes in the Cu02 planes has drastic effects on the magnetic properties, 
rortly after the discovery of high-T c cuprates, samples of La 2 _ x M x Cu0 4 with M=Ca +2 , 
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Ba +2 , or Sr +2 were studied with a wide range of probes inclu din g bulk susceptibility, NMR, 
NQR, neutron scattering. As already mentioned, neutron diffraction on powder samples 
showed that Neel state disappeared before the appearance of superconducting phase upon 
doping. 

Availability of large single superconducting crystals has made it possible to pursue the 
study of magnetic scattering into the superconducting phase. The main results of the neutron 
scattering experiments on these systems are as follows. 2D AF fluctuations are still observed 
in the doped and the superconducting samples, but with a much shorter in-plane correlation 
length, which turns out to be essentially temperature independent. As shown in Figure 1.4, 
the correlation length £ decreases from ~ 35Ato only 8Aas the Sr concentration x increases. 
The solid line, 3.8 /vTA, is the average separation between the 0~ holes in the Cu0 2 planes. 
Within rather large experimental uncertainties £ agrees with this distance quite well. The 
scattering is strongly inelastic, as in pure La 2 Cu 04 . The integrated inelastic intensity is 
comparable to that in La 2 CuC> 4 , implying that there are still sizeable moments on the Cu 
atoms in the superconducting phase [38]. 

Experiments on large single-crystals of YBa 2 Cu 3 0 6+5 in the superconducting phase by 
Shirane and coworkers [38] have revealed that magnetic fluctuations are present in these 
crystals also. The crystals were found to be orthorhombic so that the possibility of included 
oxygen-deficient tetragonal AF phases was discounted. The intensities of these fluctuations 
disappeared rapidly near room temperature. 

1.3.4 Models for the Cu0 2 planes 

As discussed above, the transport and magnetic properties suggest that these cuprate 
materials are Mott-Hubbard insulators. Very soon after the discovery of high-T c supercon- 
ductors, P. W. Anderson [39] suggested that one should start from a Mott insulator, because 
superconductivity occurs near metal-insulator transition. Moreover, he proposed that the 
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!-b&nd Hubbard model given by the Hamiltonian, 

H = —t ( a lcr a ja d" a j<r a t<x) + U (l’l) 

(ij)<r * 

ere a^(a I<r ) is the creation(annihilation) operator for spin a fermions (electrons or holes) 
site z, n ia = a\ a a tcT is the number operator, t is the nearest-neighbour (NN) hopping and 
is the on-site Coulomb repulsion, can bring out the essential physics. Also, this would 
ow us to describe superconductivity and magnetism in a unified fashion so far as the same 
ctrons are responsible for both phenomena. This suggestion is supported by the fact that 
my of the magnetic properties of the high-T c cuprates can be understood within a quantum 
iisenberg antiferromagnet (QHAF), to which the Hubbard model can be mapped in the 
•ong-coupling limit (U t). Spin- wave theory on a S = ^ QHAF predicts that due to 
antum fluctuations the ordered moment in the AF phase in 2D should be ~ 0.6/z#, which 
in close agreement with Vaknin’s measurement. The temperature variation of the spin-spin 
rrelation length in these materials agrees well with the calculated value from a 2D QHAF. 
re temperature evolution of constant energy scans in neutron-scattering measurements 
e in agreement with theoretical results obtained using dynamic structure factor given by 
rakravarty, Halperin and Nelson [40] and Tyc, Halperin and Chakravarty [41]. The energy 
pendencies of the integrated intensities of constant energy scans at different temperatures 
e also in agreement with the results obtained from the dynamic structure factor proposed 
r Chakravarty, Halperin and Nelson [40]. 

Many electron spectroscopy experiments suggest that the high-T c cuprate materials are 
.arge-transfer insulators rather than Mott-Hubbard insulators. For example, in resonant 
.tended photoelectron spectroscopy in Bi2Sr 2 CaCu 2 08 [42] the spectrum displays a high 
tergy satellite located approximately at -12.0 eV along with the main valence-band region 
tiich has a width of roughly 4 eV. The enhancement of the intensity of the satellite at 
74 eV, which is the Cu3p-»Cu3d threshold, identifies this peak with a Cu3d 9 — > Cu3d 8 
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transition [42]. Similar Cu3d 8 final state energies have been obtained in YBa 2 Cu 3 Os+$ and 
La2_ x Sr x Cu0 4 [43, 44]. In Nd 2 _ x Ce x Cu0 4 this energy is —13.3 eV [45]. Since the Cu3d 8 
peak is below the main valence-band region, this implies that the Coulomb repulsion energy 
on the Cu sites (Ud) is large compared to the charge transfer energy A between Cu3d and 
02p orbitals. This supports a CTI scenario. 

Core electron X-ray emission spectroscopy (XPS) provides additional information on the 
size of the local Coulomb correlations [46]. XPS at the Cu2p threshold has been performed 
on La 2 _ x Sr x Cu0 4 [47, 43, 48], on YBa 2 Cu 3 0 6+( s. [47, 43, 49] and on Bi2Sr 2 CaCu 2 0 8 [50]. 
Cu2p 3 / 2 core XPS shows a double-peak structure in all investigations. In La 2 - x Sr x Cu0 4 , the 
933 eV peak has been attributed to a Cu2p 3 / 2 3d 10 L final state, where Cu2p 3 / 2 denotes the 
core hole and L is a ligand hole on the surrounding oxygen sites. The satellite at 941 eV 
comprises of a Cu2p 3 / 2 3d 9 final configuration [43, 48, 47, 49, 50]. A Cu2p 3 / 2 3d 9 satellite at an 
energy larger than that of the charge-transfer state Cu2p 3 / 2 3d 10 L implies that the Coulomb 
repulsion U c d between the Cu2p 3 / 2 core hole and Cu3d 9 valence-band hole is larger than 
the charge-transfer gap. Since the on-site repulsion Ud in 3d transition metal compounds 
satisfies U d jU c d ~ 0.7 [43], one is led to a CTI picture. 

In Ols X-ray absorption spectroscopy (XAS), dipole selection rules allow only for s — )• p X) 
s — 4 Py, or s — v p z transitions depending on whether the incoming photon field E is oriented 
E\\a, E\ |6, or E | |c axis. For E\ \c a ‘prepeak’ is observed in insulating La 2 Cu0 4 and Nd 2 Cu0 4 
at E p = 530.2 eV and E p = 529.1 eV respectively. The corresponding Ols thresholds are at 
E t = 528.5 eV and E t = 528.6 eV. This feature is intrinsic to the Cu0 2 planes since it is 
absent in XAS on Nd 2 0 3 [51]. The prepeak is identified with the Cu3d 10 upper Hubbard 
band [52, 53, 54]. Since the initial state is approximately | ipi) = a \ Cu3d 9 ) + 0 | Cu3d 10 L), 
the prepeak intensity is a measure of the copper-oxygen hybridisation in the upper Hubbard 
band [52, 54]. In insulating La 2 Cu0 4 this quantity is approximately | 0 | 2 > 0.1 [52], These 
spectra are clear evidence in favour of a CTI picture for the parent compounds and cannot 
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jure 1.5: Basis (solid line bounded square) and Cu0 4 plaquette (dashed line bounded 
uare) of the Cu0 2 plane of the three-band Hubbard model. The squares and the circles 
irk the Cu 3d x 2_ y 2 and 0 positions, at energies e d and e p , respectively. The 2p x and 2p v 
bitals are located at 1,3 and 2,4, t pd denotes the transfer integral between copper and 
ygen and fj p the direct oxygen-oxygen hopping. The lattice spacing a = 2 | r |. 


' explained in terms of a LDA band structure [46]. 

Since these and other electron spectroscopy experiments (a comprehensive review of 
e topic is given in Reference [46]) support the CTI picture for the parent compounds of 
gh-T c cuprates, many authors [55, 56] propose the three-band Hubbard model, which, takes 
)th 02p and Cu3d orbitals explicitly into account, as the minimal model for describing the 
lysical properties of the Cu0 2 planes in these materials. The three-band Hubbard model 
described by the following Hamiltonian. A diagram of a portion of the Cu0 2 plane it 
presents is shown in Figure 1.5 (adapted from Reference [46]). 
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H = -tfd Y (PiA i* + h - c -) "4 E (PfarlW + h-c.) 

+£ P YPicrPi* + tdY d iAi<r 

ia icr 

+Ui £ « + C/„ 2 £ ntnf , (1.2) 

* * <*>j> 

where p{jP) and d(cfi) are the annihilation (creation) operators for holes on the O and Cu 
sites, (i,j) are nearest neighbour (NN) Cu and O sites and (l,m) are two neighbouring 0 
sites; n’s are the respective number densities. t pd is the hopping term between neighbouring 
Cu 3d and 0 2 p orbitals and t pp is the hopping matrix element between two NN O 2 p or- 
bitals. ed and e p are the on-site energies of a hole sitting on Cu 3d or O 2 p orbital. Ud and 
U v are the Coulomb repulsion terms on the Cu and O sites respectively. U v d is the repulsion 
between charges on the neighbouring Cu and O sites, a is the spin index. A = e p — e d is the 
charge-transfer gap. In the hole picture A is positive. The orbital symmetries of the type 
d x 2_ y 2 and p x>y imply additional phase factors for the hopping integrals t pd and t l J p 

t% = At pd , A = A = -1, A = A = 1, 

t% = iP’tpp, A 2 = A 4 = -1, A 2 = A 1 = 1, 

where 0 labels the central copper, and 1,2, 3, 4 the four oxygen sites on the plaquette. 

It has been argued by some researchers that the three-band Hubbard model can be 
reduced to an effective one-band model in order to describe the low-energy physics of the 
Cu0 2 planes when only the spin degrees of freedom are important. But there is no consensus, 
rather the issue is very much controversial. 


1.4 Impurity- doped cuprates 

As we have already discussed, neutron scattering experiments indicate that 2D AE 
fluctuations are present in the superconducting state, but with a much altered magnetic 
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lamics from the AF state in the parent compound. A question of fundamental importance 
io\v doping affects the electronic states of Cu +2 . Whether the electron in the Cu +2 3d x 2_ y2 
cibonding state remains localised in the highly doped samples, or it gets delocalised. The 
ifirmation of one of these possibilities was thought to be useful [57]. The particularly 
evant issues are the connection between magnetism and superconductivity, the validity of 
rious models based on localised moment, itinerant magnetism etc. 

One of the techniques used from the very early days of superconductivity to study 
ese aspects has been substitution of Cu by impurities. This ability to dope the Cu-site 
the high-T c superconductors with transition or sp elements offers an important means 
study these materials. Understanding of the effects of impurities on superconductivity 
id normal-state properties is crucial for elucidation of mechanism of superconductivity, 
re relevant aspects in substitution studies are to what degree magnetic and nonmagnetic 
ipants affect superconductivity, and how dopants with different electronic structures alter 
e normal-state transport and magnetic properties. 

Gang Xiao et al. [58, 57] of the Johns Hopkins University and Keimer et al. [4] have 
me various measurements on high-T c cuprates doped with static impurities. Xiao et al. 
ive done extensive susceptibility measurements on impurity-doped lanthanum cuprates, 
eimer et al. have done neutron scattering experiments on this system. The substitutional 
lpurities include both magnetic, such as, Fe, Ni; and nonmagnetic, such as, Zn, Al, Ga, Co; 
les. Later Mahajan et al. [2] and Walstedt et al. [59] have done various NMR experiments 
i Zn-doped Y-Ba-Cu-0 system. Impurity substitution studies in 1-2-3 systems is a little 
icky as there are two different kinds of Cu sites - in planes and in chains and one has to be 
.reful about which copper site the impurity is going into. Various observations suggest that 
.e divalent impurities like Zn +2 go into the planes, whereas, trivalent impurities, like Fe +3 , 
1+3 and Ga +3 go into the chains. There is no such complications in lanthanum cuprate as 
iere are only planar copper sites. 
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1.4.1 Effects of impurities 

All these experiments show that these impurities have four main effects on these com- 
pounds. Firstly, nonmagnetic impurities, like Zn, reduce the Neel temperature in the AF 
phase of these materials. Although they are less efficient in destroying AF order than mo- 
bile holes, nevertheless they are detrimental to AFLRO. Keimer et al. [4] measure the Neel 
temperature of La 2 _ x Sr x Cux_ y Zn y 04 single-crystals by SQUID magnetometry. They find 
that In goes down approximately linearly with increasing Zn concentration y and the T^(y) 
curve for the crystals extrapolate to zero at y -^0.26, ie., 26% of Zn doping completely kills 
AFLRO in this system. 

Secondly, of these impurities, those which substitute Cu in the Cu02 planes have 
drastic effect on T c of these materials. Measurements of Xiao et al. [58] show that in 
La 2 - x Sr x Cui_ y A y 0 4 (A=Zn, Al, Ga, Fe, Ni, Co) T c varies smoothly and monotonically with 
y. Above a critical doping level y c , T c reduces to zero. In lanthanum cuprates Fe suppresses 
T c most efficiently, with the smallest y c of 1.8%, and Ni has the largest y c of 4.2%. In con- 
trast, Zn is the most efficient in reducing T c in Y-Ba-Cu-O. In YBa2Cu 3 0 7 , replacement 
of 3% of the Cu by Zn is sufficient to reduce the T c from 90K to about 60K. 10% of Zn 
doping destroys superconductivity completely. Al and Ga are much less disruptive of super- 
conductivity in this material. This supports the site assignments stated earlier. The relative 
effectiveness of Zn in destroying superconductivity can be attributed to Zn-substitution on 
the planar sites, where it is believed that superconductivity resides. Whereas, Al and Ga, 
substituting Cu in the chains, have very little effect on T c . 

Thirdly, all these static impurities have very interesting magnetic effects in high-T c 
cuprates. In conventional superconductors magnetic impurities are most destructive of su- 
perconductivity due to their pair-breaking effect. Nonmagnetic impurities have very little 
effects on superconductivity. But, as we have seen, in YBa 2 Cu 3 07 , nonmagnetic Zn is most 
detrimental to superconductivity. Even in lanthanum cuprate, the largest y c differs from the 
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allest y c only by a factor of 2. One of the immediate questions is how the effect of dopant 
•elated to its magnetic characteristics. 

In order to study the relation between superconductivity and magnetic nature of the 
pant, Xiao et al. [58] measured the magnetic susceptibility (x) with a SQUID magne- 
neter. Interestingly, they found that irrespective of the magnetic nature of the dopant, 
T) shows a Curie- Weiss behaviour (1/T) indicative of a localised moment. All the data 
aid be well-described by the relation 


X(T) = Xo + 


3 k B (r-ey 


(1.3) 


iere xo represents the background susceptibility which does not change appreciably with 
ping, N is the doping content and p e ff is the effective moment of a dopant. In their 
alysis, Xiao et al. assumed that the local moment resides on the dopant site. As will be 
jcussed later, this is not quite correct. What is surprising about their observation is that 
the nonmagnetic dopants, that cannot carry a moment by their atomic structures, induce 
similar moment ~ 1.2ps- Fe carries a large moment of 5 p B . So, judging by their effects, 
the dopants at the Cu site are magnetic in nature. 

Nonmagnetic dopants such as Zn +2 , Ga +3 , Al +3 have inert full-shell structure (3d 10 and 
6 ). Their main role is to remove the spin of Cu +2 3d 9 (S = |). As stated already, every 
nmagnetic dopant induces a moment of 1.2 p B irrespective of its electronic structure. Such 
nmagnetic dopants do not, in general, induce moments in conventional metals. The ap- 
arance of induced moment indicates localised nature of the d state in the superconducting 
■2- x Sr x Cu0 4 . This moment, Xiao et al. argues, is nothing but the Cu +2 moment, which 
concealed from susceptibility measurement due to moment compensation in the dynamic 
r state. A nonmagnetic dopant breaks the compensation and creates a net spin of S = \ 
aerating an effective moment p e // = 9\j s ( s + 1) = 1.9 p B - The experimental p e /f value is 
>se to 1.2 pq, a 35% reduction from the theoretically expected value. This, it is believed, is 
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due to Cu3d x 2_ y 2-02p hybridisation. But, as we will shortly see, this picture of an AF plus 
a spin-hole is rather simplistic in view of the recent NMR experiments in Y-Ba-Cu-0 [2]. 

Another important observation that Xiao et al. make is that the extent of suppression 
of superconductivity has little relation with the valence of a dopant. For example, Zn +2 and 
Ni +2 , having the same valence, differ a great deal in their y c values. On the other hand, 
Zn +2 and Al +3 , Ga +3 , Co +3 having different valence states, have similar y c values. One can 
expect that the divalent impurities will not affect carrier concentration, whereas, trivalent 
ones will reduce carrier concentration with increasing doping. However, since y c is small in 
every system and insensitive to the dopant valence state, argues Xiao et al. , suppression of 
T c cannot be predominantly attributed to any change in carrier concentration. 

Recently Walstedt et al. [59] have done Cu 63 (2)-NMR studies and Mahajan et al. [2] 
have done Y 89 -NMR studies on Zn-doped Y-Ba-Cu-O. Walstedt et al. find a predominant 
T -1 temperature dependence in the Cu(2) NMR linewidth. They have explained this obser- 
vation from spin-exchange scattering of carriers from localised moments associated with Zn 
impurities in the CuC >2 planes. In their Y 89 -NMR measurements on YBa2(Cui_ y Zn y ) 3 06+<5, 
Mahajan et al. find two satellite peaks along with the main Y 89 peak. Comparing this 
with their numerical calculations, and comparing the Curie-like temperature dependence of 
the frequency shifts of the satellite lines with susceptibility data, they conclude that the 
Zn-induced local moments reside on NN Cu orbitals. This is a very significant conclusion 
as it makes the picture of the Zn-induced local moment as an AF with a spin-hole and the 
moment residing on the dopant site untenable, and calls for a more detailed microscopic 
understanding of the whole phenomenon. Since the frequency-shift of the main Y-line in 
their experiment is independent of T, they conclude that there is no change in carrier con- 
centration in Zn-doped Y-Ba-Cu-O. 

One aspect the researchers are not in consensus about is the mechanism for reduction 
in T c . Since it is seen that all the dopants, whether magnetic or nonmagnetic, give rise 
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local moments in the C 11 O 2 planes and the Fe-doped La 2 _ x Sr x Cu 04 system with the 
hest moment has the lowest y ct whereas, the Ni-doped system with the smallest moment 
j the highest y c , Xiao et al. [58] believes that an Abrikosov-Gorkov (AG) pair-breaking 
chanism is responsible for reduction in T c . Walstedt et al. [59], on the other hand, believe 
it a pair-breaking mechanism cannot explain the observed reduction in T c . Comparing 
; measured T~ l temperature dependent Cu(2) NMR linewidth with a calculation based 
RKKY-like spin-density oscillations in the Cu0 2 planes driven by exchange between the 
)bile carriers and localised moments associated with Zn-sites, they find that the reduction 
T c from the AG formula should be of the order of a tenth of a Kelvin - orders of magnitude 
taller than the observed value (~ 100K). They go on to argue that this suggests a d-wave 
.iring in the cuprate superconductors, in which case, a strong potential scatterer like Zn 
n cause pair-breaking effects. Mahajan et al. [2] from their analysis of Y-NMR linewidths 
nclude that the reduction in T c can be ~ 100K from an AG theory. They believe that 
alstedt et al. reached the opposite conclusion mainly as their y c (= C m /T), the Curie 
rm in the susceptibility in the doped systems is larger. In a recent magnetoconductance 
periment on single crystals of YBa 2 (Cui_o.oiyZno.oi y ) 30 6 +j for y — 1 and 3.5, Axnas et al. 
0] have found that phase-breaking scattering rate increases strongly with Zn concentration, 
dicating that Zn causes pair-breaking in this material. Williams et al. [61] from their Y- 
MR measurements on Y(Bai_ x La x ) 2 (Cu 1 _ y Zn y ) 4 0 8 find that magnetic pair-breaking by 
.e induced local moments can fully explain the rapid depression in T c with Zn substitution 
La-free overdoped (high hole concentration) samples. However, as the samples become 
ore underdoped, magnetic pair-breaking becomes increasingly unable to account for the 
>served decline in T c . They believe that this occurs predominantly via a reducing pairing 
teraction. 

Finally, the temperature dependence of the in-plane resistivity (p ab ) of Zn-doped cuprate 
iperconductors shows some interesting behaviour. Fukuzumi et al. [62] measured the 
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temperature variation of in YBa 2 (Cui_ y Zn y ) 3 06-f,j and La 2 _ x Sr x Cui_ y Zn y 04 over a wide 
range of Zn content and doped hole density. Interestingly, Zn substitution makes marked 
difference between underdoped and overdoped samples in these experiments. In underdoped 
YBa 2 (Cui_ y Zn y ) 3 0 6 .63 with Zn content ranging up to y 0.04 they find that is rapidly 
reduced and the compound becomes insulating after the superconductivity disappears. A 
large residual resistivity po adds to the T-dependent component. Superconductor-insulator 
(SI) transition occurs at p 0 ~ AOOyClcm. This value corresponds to the 2D resistance of 
~ 6.8AT2/Cu0 2 plane and is near the universal value hf 4e 2 ~ 6.4ATD/Cu0 2 plane. In 
La 2 _ x Sr x Cui_ y Zn y 04 for x — 0.10 and 0.15, the SI transition takes place near the universal 
2D resistance. The highly doped superconducting compounds show a contrasting behaviour. 
The x = 0.20 La 2 _ x Sr x Cux_ y Zn y 0 4 sample remains metallic even after the superconductivity 
disappears and the residual resistivity is smaller than the critical value observed for x = 
0.10 and 0.15 by a factor of 4. This shows that the underdoped compound is intrinsically 
insulating. Suppression of superconductivity drives it into an insulating state. 


1.5 Outline of the thesis 

In order to develop an understanding of the effects of nonmagnetic impurities in the 
Cu0 2 planes we study the effects of such impurities on a Mott-Hubbard AF. We find that 
each nonmagnetic impurity gives rise to a defect state within the Hubbard gap. This state 
is localised mainly on the four sites neighbouring the impurity and produces a local moment 
[63]. We also study effects of nonmagnetic impurities in a Hubbard model with only short- 
range order and find that defect states are formed in this case also. These nonmagnetic 
impurities give rise to a logarithmically diverging correction to the magnon energy in the 
long-wavelength limit in a 2D system. In an anisotropic 3D (layered) AF, this divergence 
is cut off by the interlayer hopping. We calculate the reduction in Neel temperature due to 
thermal excitation of spin waves and find that T N goes down with increasing Zn concentration 
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accordance with observations of Keimer et al. [4] in their neutron-scattering experiment 

Additional holes put in a system with static nonmagnetic impurities will go into the 
ect state as that is the lowest energy state available for a hole. So in presence of holes the 
al moment will be lost. But thermal and quantum fluctuations can transfer some spectral 
ght of the hole in the upper Hubbard band so that part of the moment can be regained. 
: present a quantitative estimate of how the moment is recovered with temperature and 
3 to quantum fluctuations. 

Defect states in a one-band Hubbard model cannot explain how local moments are 
med in systems with added holes and at low temperatures, how the observed moment 
smaller than the theoretically expected value and why all nonmagnetic impurities, irre- 
ictive of their structural details, give rise to the same moment. We attempt to explain 
;se by studying the effects of nonmagnetic impurities in a CTI described by a three-band 
ibbard model. In this case each impurity produces two defect states. The upper defect 
ite, formed slightly below the upper Hubbard band, has amplitude mainly on the four 
ighbouring Cu sites and remains occupied by electrons even in the presence of added holes 
the system. So the moment formed in this system is robust with respect to hole doping. 

We have studied the effects of magnetic impurities in a Mott-Hubbard AF. We present 
reral formalisms representing different situations. In the case impurity spin (S') equal 
the host spin (5) represented by modified hoppings around the impurity site, we find 
momentum-independent renormalisation of the long-wavelength magnon mode energies, 
the high-energy end of the spectrum, a magnon mode is split-off from the band and is 
s one localised at the impurity site. Its energy is modified because of different effective 
upling ( J') between the impurity spin and the neighbouring ones. We develop a formalism 
represent magnetic impurities by a spin-dependent potential on the impurity site as a 
tural extension of the formalism developed for nonmagnetic impurities. In this case, within 
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a T-matrix formalism, unlike the nonmagnetic impurity case, we find no defect states within 
the Hubbard gap, rather, one impurity state is formed for each spin. These impurity states 
are localised on the impurity site and the local moment in this case is because of difference 
in t- and 4,-spin densities on the impurity site. We also develop a formalism to represent 
magnetic impurities by putting multiple orbitals on the impurity site. We explicitly show 
that the Goldstone mode is preserved in presence of impurities. We also study the effects 
of these impurities on the electronic spectral properties of the system. Via a scaling study 
within a numerical HF scheme we, find a linear reduction in the magnon velocity in the long 
wavelength limit with increasing doping concentration which shows a logarithmic divergence 
in the limit of large system size. 



Chapter 2 


tatic Nonmagnetic Impurities in a 
dott-Hubbard Antiferromagnet 


The motivation for studying the effects of static impurities in a Mott-Hubbaxd AF has 
en discussed in detail in the previous chapter. In this chapter we present our formalism 
• treating nonmagnetic impurities in a Mott-Hubbard AF. The same formalism will be 
sd when we discuss nonmagnetic impurities in a magnetic system with short-range order 
hapter 3) and in a CTI (Chapter 6). 

In the first section of this chapter we briefly review the Hartree-Fock (HF) mean-field 
IF) treatment of the Hubbard model at half-filling. In the subsequent sections we introduce 
e impurity and study its effects on the electronic properties of the system. Effects on the 
agnetic properties will be discussed in Chapter 4. 

.1 HF treatment of the Hubbard model 

The quantum AF state of the Mott-Hubbard insulator has been studied in great detail 
ithin a HF-plus-fluctuations approach by Singh and Tesanovic [3]. Within HF approx- 
lation the Hubbard model on a square lattice, described by the Hamiltonian given by 
quation 1.1, has been represented in the AF state in the two-sublattice basis as, 
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-(tA e* 


( 2 . 1 ) 


where a — ±1 for f- and 4-spins, = — 2t[cos(& x a) 4-cos (fc y a)] is the free-paxticle band energy 
for a 2D square lattice and A is the gap parameter. The eigenvalues of the Hamiltonian yield 
the quasiparticle energies 


E\ = ±(A 2 + ej)>' 2 . 


( 2 . 2 ) 


The plus(minus) sign is associated with k states lying inside (outside) the magnetic Brillouin 
zone (MBZ). This shows that there are two bands for both f- and 4- spins. At the half-filled 
level, the lower Hubbard bands (LHB) for both the spins are completely occupied and the 
upper Hubbard bands (UHB) are completely empty. The system, therefore, is an insulator. 
The gap parameter is related to the sublattice magnetisation m by 2A = mU. The sublattice 
magnetisation is defined as the average over one sublattice of the expectation value of the 
z-component of the spin operator S l z = a\a z ai , where a z is the ^-component of the Pauli 
spin matrix. The quasiparticle amplitudes on the A- and B-sublattice sites are given by 
at 2 = |(1 =F aA/E t ) and bfj 2 = |(1 ± <rA/E t ). 


The sublattice magnetisation m is obtained from the difference in spin densities on a 


site: 


m = 4 E M *) 2 - ^(S) 2 ] = 4 E 


A 


Er<Ep 


N 


Er<E-p 


(A 2 + e|) x /2- 


( 2 . 3 ) 


Since 2A = mU , this leads to the self-consistency condition: 

1 - Ji V 1 
u~ (A* + el)W 


In terms of the quasiparticle amplitudes and energies, the Green’s function for the system 
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g a {k. u ) ) =< k 


\ 1 - -\h^~ ^ 

“L 

Q'kcr s^kcrs 

1 TT/T / 7*\ ^ 

to — Hq (/ c ) 

_ Q'kcrs^kas 

tin- 
kers -* 


~ E kas + i$r l 


(2.4) 


•e s refers to the two signs +/- corresponding to the upper/lower Hubbard bands. 

The spin-wave mode shows up in the form of a pole in the dynamical transverse spin 
iceptibility evaluated in the Hartree-Fock AF state. The transverse susceptibility is defined 
terms of the spin-raising and spin-lowering operators as, 


X = -?)(<& af | [S (r,t),S + (r',i') 


^af)- 


(2.5) 


ithin random phase approximation (RPA), the transverse susceptibility in the momentum 
ace in the sublattice basis representation is [3] 


x~ + (Q,ty 


x°(&fl) 

i -u x °(Q,ny 


(2.6) 


ie zeroth-order susceptibility, x°(Q, £2) is a 2 x 2 matrix in the sublattice basis and is given 
terms of the Greens function as, 


X°(Q, = &(k, u)G l {k -Q,w- n). (2.7) 

k 

ie RPA susceptibility matrix can be expanded in terms of the eigenvalues A((Q, £2) and the 
jenvectors | §a(<2, ty) of the x°(Q, £2) matrix as, 

X~ + (Q,n) = | $ A (Q, w ))<$ A (Q,f2) | . (2.8) 

mce it is clear that all the information regarding the nature of these collective spin excita- 
ms is contained in the eigensolutions {| T A (Q, £2)), A (Q, £2)}. The eigenvectors | $ A (Q, £2)) 
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represent the spin-wave amplitude in momentum space, while the spin-wave energies, Qq, 
are given by the poles at 1 — UX(Q, ft^) = 0. The wave function of the spin-wave in state Q is 
given by a superposition of plane waves on the two sublattices, the superposition amplitudes 
being given by the eigenvector $(Q,Qq), 

= (2.9) 

In the strong-coupling limit ( U » t) x°(Q, ft) has the following form up to O(^): 

x°(Q,ft) = 

where A = jj — -Jj, a = |A 2 and Bq = 

Two eigenvalues of the Xo + matrix are, 


A — Qfft 


B Q 


( 2 . 10 ) 


*/t + otCl 

with 7 q = |[cos(<5 x a) + cos(Q y a)] in 2 D. 


. _ 1 t 2 . 

A ” U A 3 1 T 


(£) 2+ i 


1/2 


( 2 . 11 ) 


where J = —■ is the NN coupling in the equivalent Heisenberg model. Equating 1 -DA max = 
0, where A max is the larger of the two eigenvalues, the spin-wave energies obtained are: 


ft Q - = ±2J(l- 7 |) 1/2 . 


( 2 . 12 ) 


2.2 Treatment of the impurity-doped system 

We now present results of our study of a Mott-Hubbard AF doped with static, nonmag- 
netic impurities. We use Hubbard model to describe 3d holes in the Cu02 planes. Zn , 
that replaces Cu^, has a filled 3d^ configuration. We model zinc-type dopants by high 
impurity potential (^ r ) off the dopant sites, to ensure the absence of any 3d hole density on 


these sites. 
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The Hamiltonian of the system in presence of impurity' is 

H ~-t {a^cLjo- + h.c.) + U ^2 n^riii + V a\ c a Ic , (2.13) 

(ij)cr i 1 i a 

ere the sum over I goes over the impurity sites. We start with the HF description of the 
doped AF host and then treat the impurity terms perturbatively within a Green’s function 
■malism. For a single impurity the perturbation can be treated exactly in terms of the 
matrix formalism. There is a caveat that this procedure is not fully self-consistent as the 
nsity modifications introduced by the high-impurity terms are not treated self-consistently. 
The Green’s function for the system with a single impurity at site I can be written as, 

= S&M + SiiHVgij + + ■■■ (2.14) 

aich can be summed up and written in terms of the T-matrix, T^w) = ^ as 

= g^u) + (2.15) 

formation regarding additional poles is contained in the T-matrix, and therefore, a study 
’ the local Green’s function in the AF state is of interest. If, for concreteness, we 

insider the impurity to be on an A-sublattice site, then we obtain, 

c f \ I s ^ (T A 

S//M -— El (2.i6) 

k k 

plot of vs. u is given in Figure 2.1 and shows logarithmic divergences at energies 

Y A 2 + (4t) 2 , a stronger divergence (w + A)~ 1//2 ln(u> + A) as ui — f —A from above, and a 
'A - In (A - uj) dependence as u -» A from below. (See Appendix A for details). The 
itersection with lfV yields the poles, and for large V we obtain an impurity state at energy 
'i = V — A (V relative to the lower band at energy -A), and a defect state inside the 
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Figure 2.1: Plot of gj f ( u>) and 1/V. 


Hubbard gap close to the upper band edge. 

The energy spectra for the two spins in presence of impurity are shown in figure 2.2 
In semiconductor terminology, such states in the gap arising from impurities are called deep 
traps , and for large V and two similar bands, occur near the mid-gap position. To reflect the 
asymptotic convergence of the trap energy to some fixed position near mid-gap as V — > oo, 
the concept of pinning has been introduced — large differences in atomic energies influence 
the trap energies only weakly [65]. However, in the AF, while the two Hubbard bands are 
identical, the two-sublattice structure is responsible for pinning the defect states to the band 
edges. 

The energy separation 77 = A - E d , of the defect state at energy E d , from the lower edge 
of the upper Hubbard band is approximately given by ^/rjlnr} ~ 1/V, so that 77 vanishes 
essentially as (1/V) 2 as V -> 00. (For details, see Appendix A). 

For the sake of completeness, we show a plot of g^(cu) and 1/V in Figure 2.3 and find 
that for f-spin there is only an impurity state and no defect state. 
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Impurity states 


N/2-1 ZZZZZZZZZZ ZZZZZZZZZZ n/2-i 

Defect state 

N/2-1 iiiiimn niinrrn 

Figure 2.2: Energy spectra for the two spins for an impurity on an A-sublattice site. 

Both the impurity state and the defect state, lying as they do outside the AF band, are 
calised. From Equation 2.15 the change in the Green’s function, 


SGij(u) = Gij(u) - gij(w), (2.17) 

le to impurity scattering can be written, for energies outside the band and very close to 
e defect state energy, E d , as: 


<>«&) »«(«,) 
SGiAw ~ Ed) = 


(2.18) 


tiere gn{u) has been Taylor-expanded near E d , and we have used 1 - Vg n (E d ) = 0 which 
.‘ter min es E d . The above equation can be rewritten in the following form, yielding the 
ifect state Greens function in terms of = gu(E d )/ \J—dg n (ui) j duj\ Ed , which gives the 
nplitudes for the defect state, 
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SG{j (iv ~ Ed) — 


- Ed + ir)' 


(2.19) 




Figure 2.3: Plot of and 1/V. 


The impurity state Green’s function can be written in an exactly parallel manner, yield- 
ing the impurity state wave function, Impurity state wave function is essentially localised 
on the impurity site in the strong-coupling limit and for large V. Continuing with the de- 
fect state properties, in the limit V oo in the strong-coupling limit, <p d has vanishing 
amplitude on all A-sublattice sites, and on the B-sublattice sites, as mentioned above, it is 
given by <p d (r) J2k 1 exp (ik.r), where r = f\ - r>. (For details see Appendix B). At large 
distances, the long-wavelength modes contribute predominantly, and this yields a power-law 
fall-off, 4> d {r) ~ (sin k c x x/x). (sin kyy /y), where k% and k% are the short-wavelength cut-offs in 
k x and k y respectively. 

The defect state density for different values of rj (or equivalently different Vs) for a 
fixed A (= 5) around the impurity site is shown in Figures 2.4, 2.5 and 2.6. We observe 
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that as 77 decreases i.e., V increases, the defect state has increasingly greater amplitudes on 
sites away from the impurity site. In other words, it is spatially more extended for larger 



Figure 2.6: Defect state density for 77 = 0.001. 


We believe that these defect states which are occupied with electrons and localised 
around the impurity sites are responsible for the presence of local moments observed in these 
systems. For a zinc atom sitting on an A-sublattice site, the defect state is occupied by an 
t-spin electron. In terms of electrons, the dopant site is occupied by two electrons (empty 
of holes), and there are (N — 1 ) t-spin and 4-spin electrons in the antiferromagnetic band 
states (in the upper Hubbard band). These electrons form compensated spin systems, and 
therefore should contribute an essentially temperature-independent (for T <C J) magnetic 
susceptibility. The lone electron in the localised defect state behaves essentially like a free 
spin in the dynamical state when there is no barrier between the up and down spin states, 
and hence should have a Curie-tvpe 1/T contribution. 


Chapter 3 


Static Nonmagnetic Impurities in a 
Spin-System with Short-Range Order 


In Chapter 2 we have discussed the effects of a static, nonmagnetic impurity in an AF 
td seen that each impurity gives rise to a defect state, which can explain the formation of 
cal moments seen in the susceptibility and NMR experiments. But moments are also formed 
compounds without AFLRO. So a study of the effects of static impurities in magnetic 
stems without long-range-order is relevant. The hope is that if we find such defect states 
this case also, then we have a more direct way of understanding the formation of local 
oments in the high-T c cuprates away from their AF phase. 

In this chapter we present the results of our study of nonmagnetic impurities in a spin- 
stem which does not necessarily have long-range-order. As in the case of an AF, we model 
nonmagnetic impurity by a high on-site potential (V) on the dopant site to ensure absence 
any 3d hole density on this site. Our main result is that defect states are formed in this 
se also. The defect states, lying outside the bands, will be localised and can give rise to 
cal moments. 

The specific model we consider is again a one-band Hubbard model. In the magneti- 
lly disordered state, characterised by short-range AF order, a completely unrestricted HF 
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decoupling of the Hubbard Hamiltonian gives [3] : 

Hhf = T- £ *![*•£]*< + £ (3.1) 

i i ^ 

where T is the hopping term, = ( a| t ajj_ ), and A, = ?7(^ r |||5 ; |^ r i). Here A t - represents 
the local magnetic field seen by an electron on site i due to electronic correlations. For the 
case of one electron per site, the (last) density term simply shifts the energy levels by U/2. 
All the Aj’s have the same norm, but they vary in direction, which is a signature of lack of 
long-range-order in the system. 


3.1 Analysis of impurities in a two-site system 

In order to gain some elementary understanding of the system, we first consider a two- 
site system as shown in Figure 3.1. On site 1 the spin is pointing along the z-direction 

(-t). 


A 

z 


1 2 



Figure 3.1: A two-site cluster. 


On site 2 the spin makes an angle 9 with the z-direction, i.e., it is tilted at an angle a — 6 i r 
from a perfect AF direction. In the atomic limit ( t = 0) the Hamiltonian [a. A] for this 
system can be written in the two-site, two-spin basis as a 4x4 matrix, 
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-A 0 0 0 

0 A 0 0 

0 0 A cos 6 -A sin# 

0 0 —A sin 9 — A cos# 


(3.2) 


Diagonalising this Hamiltonian we get two sets of degenerate eigenvalues, -A, -A and 
+A with eigenvectors 


0 
0 

V 0 / 


0 
0 

cos 9/2 
sin 9/2 J 


( 0 \ 
1 
0 

v 0 y 


and 


0 
0 

sin 6/2 
y — cos 9/2 J 


spectively. We now treat the effect of hopping as a perturbation to these states. 

We do a standard degenerate perturbation theory to first order in t. (For details see 
>pendix C). Once the hoppings are introduced, the states at —A split-off to two states at 
\ ±t cos 9/2 and the ones at +A to +A ± t cos 9/2. It also tells us that the t-spin hole 
nsity of 1 on site 1 is equally distributed to the two states at —A ± t cos <9/2. If 6 is not 
ich different from 7 r, we can assume that hoppings will transfer some spectral weight to the 
per-‘band’ which will essentially be same as that in a perfect AF system. So, in the strong- 
upling limit, distribution of t-spin weight of 1 on site 1 in the four states are 5(1 — ^), 
1 - an d respectively in -A - tcos9/2, -A + t cos 6/2, +A - t cos 9/2 and 

1 + t cos 9/2 to first order in perturbation theory. It is easy to see that to next order 
e energies in the upper and lower ‘bands’ are shifted by constant amounts, and 
spectively, independent of 9. 

From our analysis of static nonmagnetic impurities in a Mott-Hubbard AF in Chapter 2 
: know that if we place an impurity on site 1 we will have additional states within the ‘gap’ 
ien 9n (v) = 7 for some w within the ‘gap’. In the limit V — > oo this gives, 



3.2 Exact eigenstates analysis of impurities in a five-site duster 


37 


Klzjfe | KLufel i 

u> — (—A — - t cos 0/2) u — (-A — -I- f cos 0/2) 

x t 2 i t 2 

2A^ , 2A£ _ Q (3_3) 

— (A + ^ — f cos 0/2) cu — (A + ^ 4- f cos 0/2) 

For a perfect AF arrangement, i.e., 0 = x, this reduces to 


1 _ ii 

1 A 2 


W + ^ + 2A U ~ 2A 


+ 


Ji 


= 0 


(3.4) 


For oj’s within the ‘gap’, i. e. , between —A and +A, the first term is always positive and 
the second term is always negative. In the strong-coupling limit, the numerator of the first 
term is close to one and that of the second term is close to zero. So, for the second term to 
cancel the first one its denominator also has to be small. In other words, the position of the 
defect state will be close to +A. In fact, from Equation 3.4 we get, 


£j = A-fd (3.5) 

Following a similar line of reasoning one can see that for some 0 close to but not exactly 
equal to x, the defect state will be close to and slightly below the state at A + ^ — f cos 0/2 
(Figure 3.2). So for with increasing a’s the position of the defect state will be lowered. 


3.2 Exact eigenstates analysis of impurities in a five- 
site cluster 

With this elementary understanding of the effect of impurities in a spin-system which 
does not have a perfect AF order, we proceed to present our analysis of nonmagnetic im- 
purities in a five site cluster of Hubbard model, within an exact eigenstates analysis. The 
system we consider is shown in Figure 3.3. 
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+A+t cos(0/2) 

t ^ 

+A- 1 cos(0/2) + ~a 


+A 


-A+ 1 cos(0/2) — 4r 

2A 

-A- 1 cos(0/2) — — 
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Figure 3.2: Position of the defect state in a two-site cluster. 



Figure 3.3: A five-site cluster. 
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We assume that on the central site the spin is pointing ‘up’. On the four neighbouring 
sites the spins are tilted at an angle oc with respect to a perfect AF arrangement, i.e.. 
they make angles 9 — ir — a with the z-direction. In the magnetically disordered state, 
characterised by short-range order, that we are talking about, a is a measure of the inverse 
correlation length of the system, in fact, a ~ 7t£ - 1 [66]. We now find the position of the 
defect state as a function of a or the inverse correlation length. 

In the site-plus-two-spin basis, the Hamiltonian matrix for the system is a 10x10 matrix 
as shown in Appendix D. 

The energy-scale in the problem is set by taking t — 1. We take A = 4.0 and V = 15. 
Again, if a nonmagnetic impurity is put on site 3, the position of the defect state will be given 
by that to for which G 33 {uj) = 1/V. So, the task is to calculate G 33 {u) from the eigenvalues 
and eigenvectors of H for some chosen 9 and to find that value of to for which G 33 (li) = 1/V. 
By the construction of [. H ] , the f-spin amplitude on the central site for the n-th energy state 
is given by ^ n (5), where [-0] *s are the eigenvectors of [H]. That for the 4.-spin is given by 
-d n (6). So, the local Green’s functions on site 3 are given by, 


— IZ 


£33 M — ZZ 


^n(5) 2 

A„ 

•0n(6) 2 


~n ^ ^ n 


(3.6) 

(3.7) 


A n being the n-th eigenvalue of the Hamiltonian. It turns out that for u > within the ‘gap’, 
the 4,-spin Green’s function is negative so that ^33(0;) can never be equal to 1/V for positive 
V. We solve Equation 3.6 for oj to see if there are any states in the gap. 

We find that with increasing a the five states in the upper ‘band’ spreads out in energy. 
But there is always a value of ui below the lowest of them which satisfies Equation 3.6. So 
irrespective of the range of order in the system, we always have a defect state within the 
‘gap’. Figure 3.4 shows 77, the position of the defect state from A as function of inverse 
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•relation length a of the system. 



Figure 3.4: Position of the defect state with changing a. 


The reason why 77 increases with a can be understood as follows. As we have discussed 
. Chapter 2 , the two sublattice structure with unequal distribution of the two kinds of spin 
ansities on them is responsible for the defect state to be formed near the band edge. With 
icreasing a the difference in densities of the two spin species will be lowered and hence the 
-feet state will be formed more and more near the mid-gap position. 





Chapter 4 


Impurity-Scattering of Spin- Waves in 
a Doped Mott-Hubbard 
Antiferromagnet 


In this chapter we extend our T-matrix analysis to examine impurity scattering of spin 
waves within the Hubbard-model representation of an antiferromagnetic insulator doped 
with static, non-magnetic impurities. The transverse spin-fluctuation propagator has been 
obtained within the Random Phase Approximation (RPA), and impurity-scattering-induced 
correction to spin-wave energy of long-wavelength modes has been obtained to first order 
in doping concentration x. In agreement with earlier results [5], we find a logarithmically 
divergent reduction in the spin- wave energy of long- wavelength modes in two dimensions. 
We also extend this analysis to a layered Mott-Hubbard system and show that the impurity- 
scattering-induced softening of spin-wave energies has a direct bearing on the almost linear 
reduction of the Neel temperature with x as seen from neutron-scattering studies [4]. 

At the RPA level the transverse spin-fluctuation propagator for the doped system is given 
by [x -+ (fl)] = [x 0 (O)]/l-Cf[x°(n)], where [x°(0)] is the bare, antiparallel-spin particle-hole 
propagator, and the pole 1 — £/[x°(fi)] = O' yields the spin-wave energies. We have obtained 
the correction to spin- wave energies perturbatively in terms of [dy 0 ] = [y°] - [XhostL the 
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rturbation introduced due to impurity terms in the host antiferromagnet. \5x ] has been 
duated, in turn, in terms of the modification 6G a to the one-particle Green’s function 
tained within the T-matrix analysis. 


.1 Spin-wave spectral properties of a doped 2D an- . 
t iferr omagnet 

In this section we focus on the spin-wave spectral properties of a doped, two- 
mensional AF. In order to obtain the impurity-induced perturbation [5x°(G)]y = 

"fT - Q) - — fl)| , mentioned earlier, we consider first the lo- 

.1 density modifications caused by the high impurity potential. 

For simplicity we consider the host antiferromagnet in the strong coupling limit (U t). 
i this limit, to order (f 2 / A 2 ), the hole and particle densities on sites of the two sublattices 
and B have values 1 — f 2 /A 2 and t 2 /A 2 . Suppose on A sublattice sites, the spin-f and 
>in-i hole densities (coming from the lower Hubbard band) have values 1 — t 2 / A 2 and t 2 / A 2 
spectively. If now on a site I € A, an impurity is placed with a high impurity potential 
' -4 oo), the small local spin-4, hole density (t 2 / A 2 ) is displaced, mainly to the four nearest 
hghbour (NN) sites. This can also be seen by explicitly evaluating the correction to spin-4, 

3le density on the NN sites from the Green’s function in Equation 2.17. If J labels the NN 
tes and 8 = fj - fy, then for 8n l j = —i f %%8 Gjj(uj), we obtain, using Equation 2.17 and 
Durier transforming the host Green’s functions, 



he two host Green’s function terms in Equation 2.17 contribute the two e £ / 2A factors, being 
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off-diagonal in the two-sublattice basis. In the V — > co limit the denominator 1 — V 
has been approximated by —Vgfj(u), as the host Green’s function is nonvanishing for u > 
lying within the bands. Also we have used Ylk e k e k-Q =7 q = (cos Q x + cos Q y )/2. Here the 
frequency integral has been performed by replacing with a contour integral over the upper 
half plane so as to pick poles which correspond to holes in the lower Hubbard band, and 
only terms up. to order t 2 / A 2 have been retained. This correction to spin-4 hole density on 
NN sites can be equivalently considered in terms of a correction to the local spin-4- spectral 
function in the lower Hubbard band. In the same manner we obtain the correction to the 
local spin-t spectral function in the upper Hubbard band, yielding an identical correction to 
the spin-t particle density on NN sites, again due to displacement from the impurity site. We 
are now in a position to evaluate the correction [<5x°(^)]ij = [x°(^)]ij — [Xhost(^)]»j> to the 
bare, antiparallel-spin particle-hole propagator. From the correction to one-particle Green’s 
function given in Equation 2.17, we obtain: 

[ix°(Sl)]« = « / fl [4(aO«J}i(w - a) + iG&Ms&w - £2) + SGlMSG^w - n)] (4.2) 

Now, in the V — * oo limit, all extended-state wave functions have a vanishing amplitude on 
the impurity site. The site-localised impurity states are at infinite energy and therefore do 
not contribute to the particle-hole process. Therefore [x°(^)]tj must vanish if i/j = I, i.e. 
the Nth row and the j- th column of the [x°] matrix are completely zero. This directly leads 
to the result given below, which can also be obtained from Equation 7.4 by substituting 
SG^ = —gfj for i/j = I, which similarly follows from = 0 for i/j — /, 


Px°(n)fe = -IxLMh (i/j = *) ( 4 - 3 ) 

For i/j — I, we have thus obtained the correction [5x°(^)]q in terms of the correspond- 
ing matrix elements for the host antiferromagnet. The undoped (host) antiferromagnet 
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s been studied eaxlier in detail, in both the strong and intermediate coupling regimes 
67]. Here we consider the strong coupling limit which is the simplest to deal with ana- 
tically. Translational invariance holds only within the two-sublattice basis, and therefore 
Lt(^)k = (l/A r )EQ[Xhost(Q^)]e i<5 - (rl_fri) - where [\i ost (Q, ft)] is a 2 x 2 matrix in the 
/o-sublattice basis, and is given by Equation 2.10. 

Lastly, we consider the diagonal matrix elements on the NN sites. From our 

irlier discussion of the displacement of particle and hole densities by the strong impurity 
otential, a straightforward calculation utilizing the corrections to spectral functions yields 
X°(f l)]jj = f 2 /4A 3 . The various matrix elements of the correction [<5\ ,0 (fi)] to the bare, 
atiparallel-spin particle-hole propagator are summarised below: 

[■5x°(f2)]„ =-^ + |j(l + 57 ) 

[«x°(n)]» =j|r 

(4-4) 

We now turn to transverse spin-fluctuation modes in the doped system and study the 
npurity-scattering induced renormalisation of spin-wave energy. As mentioned earlier, we 
btain corrections to spin- wave energy perturbatively by treating [Jx 0 (^)] = [x°(^)j — 
thost(^)] as ^e perturbation matrix, and determining corrections to eigenvalues of [Xhost(^)]- 
’he eigenvalues and eigenvectors of [Xhost(^)] can he labeled by wave- vector Q, and obt ain ed 
i terms of the eigensolutions of the 2 x 2 matrix [x£ 0St (Q, ft)]- If IQ) = (a# (3 q) represents 
he eigenvector of [Xhost( < 3> ^)]> so that ctQ and ,6q are spin- wave amplitudes on A- and B- 
ublattice sites respectively, then the eigensolutions {A^, \Q)} of [Xh 0S t(^)]> to second order 
i $1/2 J) and Q, are given by: 


A (oj _ i t 2 ( Q 2 



2 


(4.5) 
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Before proceeding with the perturbation theory we introduce a transformation which 
renders the perturbation theory convergent, and therefore easier to implement. This trans- 
formation pertains to the freedom in choosing the local matrix element [x°(f 2 )]/j arbitrarily. 
Since the impurity site I is completely decoupled from all other sites in the V -> oo limit, 
eigenvalues of [x°(^)] having extended eigenvectors would be independent of the impurity- 
site matrix element [x°(fl)]//- Therefore, except for the impurity-site localised eigenvector, all 
extended-state eigensolutions will remain unchanged with reselect to variation of [x°(^)]i7- 
We choose to shift [x°(fi)]/j up by 1/17, so that [<5x°(fl)] vanishes to order 1/U , which results 
in a fast convergence of the perturbation theory and, as we shall see, to determine corrections 
to the spin-wave energy in the long-wavelength limit, we do not need to go beyond second 
order. We now proceed with the perturbative evaluation of correction 5A<g(fi) to eigenvalues 
of [Xhost(^)] due to perturbation [$x°(fi)]. From perturbation theory we have, 


<5A Q (fi) = (Q|5x°(^)IQ) + E 

Q' 


|<Q|£x°(tt)IQ')l 2 


X ( 0 ) 

a Q' 


Now, for a single impurity on site I (in A sublattice) the required matrix elements are 
evaluated using [$x°(fl)] from Equation 2.10, and the spin-wave amplitudes atg,0Q given in 
Equation. 4.6. Summing over all NN sites of I and using 7 q = we obtain after 

substituting for the spin- wave amplitudes, and evaluating the matrix element to first order 

in Q and Q : 

mx°\Q‘) = 

(1 + A) aga# + A 5 % 7 < 3 _ 0 - + «q/%.7q' + 

= — — [—1 

N A 3 L2JJ 


(4.8) 
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nee we axe only interested in 0(ft 2 , Q 2 ) corrections in 5Xq(CI) in order to obtain the spin- 
ive velocity renormalisation, we can ignore higher-order terms in (<2|5x°l<2 / ) in second 
der correction dAg^(I2). For the same reason higher-order terms in the perturbation series 
in be dropped if we are only interested in the low-energy and long-wavelength modes. For 
le first-order correction SAq } (Q), we do need to consider {Q\5x q \Q) up to second order in 
2 and Q 2 . We find that the order Cl 2 terms cancel exactly, so that for a single impurity on 
site I of the A sublattice we obtain 


(Q\5x°mQ) 


1_ f_ 

N A 1 


2J 2 


(4.9) 


Now we extend these results for a finite impurity concentration x A /x B of impurities on 
ae A/B sublattice. For an impurity in the B sublattice Cl is replaced by —Cl. If x A and z B 
re the impurity concentrations on A and B sublattices respectively, then summing over all 
npurity sites I, we obtain for the first and second-order corrections: 



Cl fx A ~x B \ (Q 2 \ 

2j{ 2 J \2 ) 


^ Xa +Xjb j 



E E 

ieA/B I'eA/B 


NO 04) 

A°) >,(0) 

a Q - a Q‘ 


e -z(Q-0').(r/-r>) 


(4.10) 


(4.11) 


If we now perform configuration averaging with = x A /x B &iv (since 

? ^ Q'), so that terms in the product (Q|dx°|Q')(Q'|dx°|Q) pair-up with the same impurity 
ite, then the configuration-averaged second-order correction (dA^ 5 ) is obtained as: 



(4.12) 


a the symmetric case when x A = x B = x, the 0(CI) term in (dA^) drops out, and substitut- 
ag for A^ — Aq) from Equation 4.5 in Equation 4.12, we obtain the following expressions for 
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the first- and second-order corrections, and therefore, up to second order, for the eigenvalue 
A q(O) we have: 



The spin-wave energy Qq can now be obtained from the pole 1 — UXq(Q.q) = 0. To first 
order in impurity concentration x , this yields for the spin- wave energy CIq of long- wavelength 
modes, 


Q,q = 


V2JQ 
= V2JQ 


1 ~ x ~ 4x T7 E 


N^Q' 2 -Q 2 \ 


1-x-Ax 

U Q 2. 


(4.16) 


where Q c is an upper momentum cutoff of order 1. Thus in two dimensions there is a 
logarithmic singularity as Q -> 0, indicating a logarithmically divergent reduction to the 
spin-wave velocity. Impurity scattering also introduces an imaginary part in the spin-wave 
energy leading to a damping of the spin-wave modes. 


4.2 Spin-wave spectral properties of a doped layered 
antiferromagnet 

We now extend our analysis of impurity scattering of spin waves to a layered anti- 
ferromagnet. All high-Tc cuprates show 3D antiferromagnetic long-range order (AFLRO) 
below the Neel temperature, due partly to the very small effective interlayer coupling be- 


s 
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zeen copper-oxide planes arising from the orthorhombic distortion. The small in-plane and 
lt-of-plane anisotropies which lead to anisotropy gaps in the spin-wave spectrum also con- 
ibute to the finite Neel temperature [68]. In doped cuprates wherein impurities like Zn 
:place the planar Cu sites, T N is found to decrease almost linearly with increasing impurity 
oncentration [4]. We have therefore examined effects of impurity scattering of spin waves on 
lagnetic dynamics due to thermal excitation of spin waves in the doped system, aiming at 
quantitative understanding of reduction of Neel temperature with impurity concentration. 

If r denotes the ratio of the interplanar and planar hopping strengths, then the free- 
article energy for the layered system gets modified to = — 2f[cos k x + cos k y + r cos k z ], 
.gain retaining terms up to order (f 2 /A 3 ) in the bare antiparallel-spin particle-hole propa- 
ator, we have [3] 


[xLt(4fi)] = ^l-^3(l + r 2 /2) 


L7q 


7Q 

1_H 

2 J 


(4.17) 


’'here the spin- wave energy scale J is related to the planar exchange energy J = 4t 2 /U by 
= J(1 + r 2 / 2), and = (cos Q x + cos Q y -I- r 2 cos Q z )/{ 2 + r 2 ). The spin- wave energy has 

he same form Qq = 2 J — 7 ^ . For high-T c cuprates r <1, so that for long- wavelength 
lanar modes we have 7 J ~ 1 ~ Q 2 p /2 - r 2 ( 1 - cos Q z ), where Q p = ^Q 2 + Q 2 is the planar 
lomentum. The zeroth-order eigenvalue in Equation 4.5 therefore gets modified to: 


< n > = V~ ^ 1+r2 M (f - 5 (§) 2 + y(l - cosQj) (4.18) 

nd up to second order in Q./2J , Q p and r the spin-wave amplitudes remain unchanged as 
iven in Equation 4.6. The impurity-induced perturbation [$x°(J2)] to the antiparallel-spin, 
article-hole propagator can be evaluated using the analysis presented before, and the matrix 
lements given in Equation 4.4 are modified to: 
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(TO, =-I + jL(i +r v 2) ( 1+ !) 
fo°(n)]/j = [i/(n)] w = ~ 

[<5x°(^)]/k = [ 6 y Q ( 9 .)} KK = ) A r - (4.19) 

where J and K are respectively the in-plane and out-of-plane NN sites of I. As expected, the 
out-of-plane couplings are reduced by a factor of r 2 . Perturbative corrections to eigenvalues 
of [Xhost(^)] can now be calculated as before using these matrix elements of the perturbation 
matrix [5y°(f2)]. For a single impurity on the A sublattice, in terms of the unchanged spin- 
wave amplitudes (»q (3q) given in Equation 4.6, we obtain for the matrix element between 
spin- wave states |Q) and \Q'), 

(Q\6x°\Q') = ~(l+r 2 /2)x 

^1 + — y ^ (XqOCqi 4- PqPq'Iq-Q' + 0c qPq i 'Yq i d" 0Q a Q'7Q e Q) I ^ (4.20) 

which has the same form as Equation 4.8, however, in terms of 7 q defined earlier. For 
Q' ^ Q we need to evaluate this matrix element to first order only, and substituting the 
spin-wave amplitudes we obtain the same result as in Equation 4.8 except for the additional 
factor of (1 + r 2 /2). For the first order correction we do need to evaluate the matrix element 
(Q\$X°\Q) to second order in Q, Q p , r as before, and we obtain 

<«I*X°IQ> = 4iL( 1 + r 2 /2) [A + f + r 2 (l - cos Q,)] . (4.21) 

Now for a finite impurity concentration x in the symmetric case (x A = x B = x), with 
impurities randomly distributed on sites of the two sublattices, summing over all impurity 
sites as before and performing the configuration average, we obtain for the first- and second- 


order corrections: 



) 
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(<> = £(1 + r V 2 ) 


r o ! 


R + r 2 ( 1 - cosQj) 


x. 


(4.22) 


Q\ 2 1 


A 3 ( 1 + r / 2 )( 2 j) n (Q (, 2 — Qp) + 2 r 2 (cos Q z — cos Q' 2 ) 


- x. (4.23) 


'herefore up to second-order in perturbation theory, where terms up to 0 (fi 2 ,Q 2 ,r 2 ) have 
een retained, the corrected eigenvalue is given by: 


A e( n >=(J-|j( 1 + r2 / 2 >* 

+ y(l — cos Q z )^ (1 — 2x) — ^ ( 57 ) (1 “b 8 a: S r (Q p , Q z )) 


(4.24) 


here S r ((3p,Q,) denotes the sum ( 1 /JV) [(Qj , 2 - Q 2 ) 4- 2 r 2 (cosQ z - cos Q',)] -1 . The 

)in-wave energy Qq is now given by the solution of 1 — UXq(Q)=0, and in terms of the 
)in-wave energy Q, 0 q = \/ 2 J[Q 2 + 2 r 2 (l— cos Q Z )} 1/2 for the undoped layered antiferromagnet, 
e obtain: 


Q 0 = Q,° q [1 - x - Ax S r (Q p , Q g ) ]. (4.25) 

We now examine the impurity-scattering-induced renormalisation of spin-wave energy 
)r the layered system. We show that the interlayer hopping term acts to cutoff the loga- 
thmic divergence observed earlier in two dimensions. For this purpose we study the sum 
r(Q p , Q z ) in different limits. Replacing the sum over Q' by integrals over Q' z and Q p , and 
lcluding a factor of 2 to account for the identical contribution from the corners of the Q p 
rillouin zone, we obtain: 


Sr(Q P ,Q z ) 


2 r 191 f Qc 1 

i-n- 27r Jo 2tt (Q' p 2 -Qj) + 2 r 2 (cos Q z - cos Q' z ) 
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i_ Qt 

2tt J-* 2t r \2r 2 (cos<3- - cosQ'J - Q 2 


(4-26) 


where Q c is a planar-momentum cutoff of order 1. Since we are interested in the temperature 
regime kT <C J and predominantly long-wavelength spin-wave modes with energy JQ P < kT 
are thermally excited, therefore we are really interested only in long-wavelength modes with 
Q P 1- Consequently the cutoff Q c is typically much larger than Q p and r, and therefore 
an additive 1 in the argument of the logarithm in the above equation has been dropped. We 
now examine S(Q P , Q z ) in two limiting cases with respect to the planar momentum Q p and 
the hopping ratio r. 

(a) Q p r: In this limit the r 2 term can be dropped and the integral simply yields 
S r (Q p , Q z ) rs (l/27r)[ln(Q 2 /<3 2 ) + wr], so that the planar momentum itself acts as the cutoff, 
just as for the planar system, and the spin- wave energy reduction goes as In {Q C [Q P ). 

(b) Q p -C r: In this limit the Q p term can be dropped in the integral and in the limit Q z —¥ 0, 
we obtain S r (Q p , Q z — ¥ 0) = (1/27 r) ln(0 2 /r 2 ). The logarithmic divergence observed for the 
two dimensional system as the planar momentum goes to zero is therefore removed by the 
interlayer hopping term as expected, and the spin- wave energy reduction goes as ln(Q c /r) 
in this limit. We can therefore write the renormalized spin- wave energy as, 



We have neglected the non-singular correction, and a is a constant of order unity which has 
a value A/tt in the two limits considered. 
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.3 Reduction in Neel temperature 


We now turn to the magnetic dynamics in a doped layered antiferromagnet and examine 
ects of impurity scattering of spin waves on reduction in sublattice magnetisation due to 
ermal excitation of spin waves. Qualitatively, the softening of spin waves due to impurity 
attering causes an enhancement in their thermal excitation at finite temperature, which 
ads to a greater thermal reduction in sublattice magnetisation, and hence a lowering of the 
del temperature. Using the expression for the renormalized spin- wave energy, we obtain the 
duction in sublattice magnetisation with temperature within the linear spin-wave theory, 
xtrapolation of this to the temperature where the sublattice magnetisation vanishes yields 
quantitative estimate of the Neel temperature. This estimate ignores the temperature 
gion where critical fluctuations are important. From spin-wave theory the reduction in 
iblattice magnetisation at a temperature T is given by, 


•6M(T) = 


1^2 J 2 

A r _ Tin e^Q -1 
Q v 

dQ z rQ*Q p dQ p 2J 


= 2 r&r 

J- 7T 27T Jo 


2 7T Qq - 1 ’ 


( 4 . 28 ) 


here the spin-wave energy can be written as: Qq — Uq[1 - axln(Q c /(Q p , r) max )]. Now, 
the spin-wave-energy expression the impurity scattering term is ln(Q c /r) for Q p < r. 
owever, we can approximate this term as ln(Q c /r) in the whole range of the Q p integral, 
his is because thermal excitation of short-wavelength spin-wave modes ( Q p rsj Qc ~ l) 
suppressed at low temperatures ( kT <C J). Therefore only long wavelength modes are 
-edominantly excited, so that Q p is closer to r than to 1. Hence the weak logarithmic 
jpendence over Q p can be neglected in the first approximation. Thus within the q p integral, 
3 can approximate the spin-wave energy expression as Qq = (^[l-cur ln(Q c /r)], so that the 
lpurity-scattering-induced renormalisation is seen effectively as a momentum independent 
normalisation of the spin- wave energy scale. 
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The finite-temperature reduction in sublattice magnetisation arising from thermal ex- 
citation of spin waves in the undoped layered antiferromagnet has been studied earlier [69], 
and the integrals over Q p and Q z have been evaluated in various temperature regimes. For 
kT 2 Jr, which is about 20K for La 2 Cu0 4 [68], an essentially linear falloff of sublattice 
magnetisation with temperature was obtained: - 5M°(T ) ~ (kT/J) ln(kT/Jr). Incorpo- 
rating the renormalisation due to impurity scattering, and neglecting the weak logarithmic 
temperature dependence, we obtain to first order in x: 

-5M(T) (l + ax\n-) (4.29) 

J r \ rj 

Now a quantitative estimate for the Neel temperature can be obtained by extrapolating this 
linear fall-off to the temperature axis which yields the temperature at which the sublat- 
tice magnetisation vanishes. Taking the zero-temperature magnetisation as 1.0 (neglecting 
quantum spin fluctuations which reduce M( 0) from 1 to nearly 0.6), we obtain the Neel 
temperature for the doped antiferromagnet by setting — SM(T N ) = 1, 

T N (x) ~ T n ( 0) (l - ax In ^ (4.30) 

where Tn( 0) = J/ln(l/r) provides a quantitative estimate of the Neel temperature for 
the undoped layered antiferromagnet. In La 2 Cu0 4 the ratio r of interlayer and planar 
hopping terms is ~ 6 x 10 -3 [68] so that the ln(l/r) term is about 5. This indicates, upon 
extrapolation, that the Neel temperature vanishes at impurity concentration of roughly 20%, 
which agrees well with experiments. It must be mentioned that in La 2 Cu0 4 there are in- 
plane and out-of-plane anisotropies as well which lead to anisotropy gaps in the spin-wave 
spectrum. Anisotropy gap by itself also provides a cutoff in the logarithmic singularity 
in two dimensions. Therefore a more realistic study for La 2 Cu0 4 must incorporate both 
the interlayer hopping and the anisotropy gap A«. For the undoped system the magnetic 
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mamics involving both these has been studied [68], and roughly \Jr 2 4- (A a /2J) 2 is seen to 
place r as the effective cutoff. Since in La2Cu0 4 both r and A a /2J are of the same order 
j 6x 10~ 3 ), the falloff of with x will be nearly same as in Equation 4.30. 




Chapter 5 


Local Moment in Presence of Added 
Holes 


We now consider doping additional holes into an impurity-doped Mott-Hubbard AF. 
We consider the simplest case where one hole is added in presence of a single impurity. The 
defect state being the lowest in energy for the holes, the hole will go into this localised state. 
So in presence of the hole the local moment will be lost. But fluctuations - both thermal 
and quantum - can transfer some of the spectral weight of the trapped hole to the UHB. 
In other words, now there will be some spectral weight of the electrons in the defect state, 
'ensuring that part of the moment is recovered. In this Chapter we present an estimate of 
how the spectral weight of a trapped hole is transfered to UHB and at least part of the local 
moment survives even in presence of added holes. 

In order to calculate the recovery of moment due to thermal fluctuations, we use a 
numerical Hartree-Fock calculation on a finite-sized lattice. We begin the Chapter by giving 
a brief outline of our numerical studies on a finite sized Hubbard model in the first section. 
The later sections give estimates of the recovered moment. 
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1 Numerical HF studies on a finite-sized Hubbard 
model 

In this section we present the general scheme for numerical HF studies on a finite-sized 
lbbard model. The model is again described by the Hamiltonian 

H = - ( a L a jv + h - c -) + Uj2 n i t n 4> ( 5 - 1 ) 

(b').o' * 

,d the technique we use is the same as in Reference [70]. 

In the HF appro xim ation relevant to symmetry breaking along the z-direction only, the 
teraction terms can be written as, U (n it )n,- 4 . + U (n iir )n^, so that the on-site energy of spin 
= — a hole on the i-ih site is U{n ia ). In the site basis, matrix elements of the Hamiltonian 
r spin a are given by, (i | H a | i) — U (n l5 -) and (i \ H a \ j) = —t where i and j are NN 
;es. Choice of an initial configuration of spin a densities on the lattice gives the starting 
amiltonian matrix H c which is diagonalised to give the eigensolutions E ia , 4>i a . From these 
,e spin densities are evaluated using (n ia ) = ^s Jcr <£ F (di a .) 2 . Now H d is formed from these 
:nsities and diagonalised, which, in turn, yields (n^). The procedure is continued till 
lf-consistency is achieved. The impurity is modeled, as before, by putting a high on-site 
jtential V on the impurity site for both the spins. 

In our model the energy scale is set by t, and we take Ujt = 8 and V/t = 15. We take 
10 x 10 system with periodic boundary conditions. In the pure case there are two bands, 
ich containing 50 states, for both f- and !-spins. For each impurity there is one hole less 
impared to the half-filled case. So with a single impurity on an A-sublattice site there are 
) t- and 50 |-spin holes. HF calculation in this case gives rise to a defect state and an 
ipurity state for the f-spin and only an impurity state for the 4--spin just as we saw in 
hapter 2. Now there are 49 states in the LHB for the t-spins and the 50th state is the 
5fect state. This defect state, since it is occupied by electrons, we have argued, can give 
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rise to local moment. 


5.2 Recovery of moment due to thermal fluctuations 





Figure 5.1: Occupation number of electrons in the defect state. 

If an t-spin hole is now added, it will go in the defect state. In this section we study 
how the spectral weight of the hole is transfered to the UHB with increasing temperature. 
We consider a 10 x 10 Hubbard model with a single impurity as discussed in the previous 
section. The energy difference between the defect state and the UHB, denoted by 77 , is found 
to be 0.189f for the set of parameters chosen. At zero temperature, the density of electrons 
in the defect state is zero. To calculate the density of electrons in this state with increasing 
temperature, we find out the Fermi energy, n, of the system at different temperatures. This 
is done by equating £i where e i’ s are the ener §y eigenvalues of the t-spins, to the 

total number of t-spin holes in the system, he., 50. The Fermi factor 1+e (E d - fl )/H- then gives 
the occupation number of holes in the defect state. So the number of electrons in the defect 
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:e is simply 1 — ^ > which is also a measure of the local moment, the variation 

which is shown in Figure 5.1. Since the only temperature scale in the system, as far as 
:rmal excitation of holes to the UHB is concerned, is set by tj, we plot the occupation 
iiber as a function of As expected, the occupation number of the electrons in the 
bet state is small at low temperatures as a small fraction of the defect state hole is excited 
the UHB. Then it grows with increasing temperature and saturates at higher It is 
lieved that the one-band Hubbard model in the intermediate-coupling (U ~ free-particle 
ndwidth) regime can represent the CuC >2 planes in the high-T c cuprates. If the parameters 
: have taken are realistic, (f = 1 eV and U = 8 eV) then a temperature of 100K corresponds 
££ = 0.04 and we find that even at this temperature the local moment present in the 
feet state is quite small. 


.3 Recovery of moment due to quantum fluctuations 

Due to interaction of the holes with the symmetry-restoring spin waves, part of the 
feet state spectral-weight will be transfered to the UHB. So, beyond mean-field, at the 
in-fluctuation level, part of the local moment can be recovered in presence of added holes 
en at zero temperature. For analytical simplicity we consider the case when impurity 
>tential V is not too large so that the defect state wave function is essentially localised 
i the four NN’s of the impurity site. For a single impurity on an A-sublattice site, we are 
terested in processes like: 

X" + (O) 

✓ s 

/ V 

/ \ 
t \ 

/ V 

I \ 

I I 

I I 

> 1 J* 2 

<ta-E d ) j ; Gl(G> " 0) * df(co- E d ) 
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where G t (E d )'s are the local Green’s functions on the NN sites to the im purity, J x and J 2 are 
two NN s of the impurity site residing on the B-sublattice. The self-energy correction due to 
such interaction of the quasiparticles with the spin waves and the resulting transfer of the 
spectral-weight to the UHB has been calculated in detail for a pure AF [3]. We proceed along 
similar lines here. In the sublattice basis and momentum space, the self-energy correction 
due to the process is 


e(*,w) = ^£ / |h^V + ( 0 ,n)G*(£-«,w-n), ( 5 . 2 ) 

Q 

where it is understood that the corresponding elements of the two matrices, x and axe to 
be multiplied. We note that G l (k—Q, ui—Cl) has to be in the upper band as all the lower band 
4,-spin states are occupied. The x~ + matrix has the following form in the strong-coupling 
limit [3]: 


X~ + (Q,n) = - 


1 

«ia 

1 

r~H 

l 

_ 7 _ 

1 

1 

1 

2 

. -Tq 

1 

ala 

+ 

r-H 

V 1 “ 

H — Qg + it] 

Q + Qq — it]_ 


(5.3) 


where O^’s are the spin-wave energies. Now picking up the appropriate matrix elements, 
replacing the G-integral by a contour integral and taking residues at the spin-wave poles we 
have, 


Ebb(^, u;) 



f 2 /A 2 


2J ' V 1 “ 7 ! w 


E U+ i5 ' 


(5.4) 


As the 4,-spin Green’s function is in the upper band, only one of the spin wave poles, which 
has a negative imaginary part, contributes in the contour integral. It has been assumed that 
the spin wave amplitude on the sites neighbouring the impurity remain essentially unchanged 
in presence of impurities. So for two particular NN sites J\ and J 2 , 
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EjWaM = 


U 2 1 




2 N . 

fc Q 


0$ 

2J 


t 2 / A 2 




w ^<3 -^f-Q + ^ 


^k.{fj x -r-j 2 ) _ ( 5 _ 5 ) 


In the strong-coupling limit taking E®_^ ~ A, so that there is no /c-dependence, and 
ing the fc-sum yields, 


- 2 A' ? 

Q 


1 + 


fV\ 


t 2 /A 2 


2 J 




a; 


fig — A + 2(5 


- fj 2 ). 


(5.6) 


lis shows that the self-energy due to interaction of the quasiparticles with spin-waves is 
agonal in position space. 

The defect state weight of 1 is equally distributed on the 4 NN sites of the impurity. So, 
corporating the self-energy correction due to interaction with spin waves, the local defect 
ate Green’s function, denoted by Q is modified by 


5G]j[u ~ E d ) = - — A_s( w ) A- (5.7) 

uj — Jb d — td uj - Ed — i6 v ’ 

dding contributions from all the four nearest neighbours, the portion of the spectral-weight 
ansfered to the UHB is 


4 X Le„ ~ E *)- (5.8) 

straightforward integration gives this to be Ta(i + ^ 1 

2 v N y 1 2J J 7^2 ,2 * 

V ^ 

This is correct only in the limit of small V , when the assumptions that 77 » E anc j that 
16 defect state “ Realised only on the four NN’s of the impurity are valid (see Appendix B). 
Iso, the first-order perturbative calculation we have done is meaningful only in the limit 
small - 5 . At this stage we cannot perform the Q-sum as we do not have an analytical 
cpression for the spin-wave energy in presence of impurities over the entire Q-range. 
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Static Nonmagnetic Impurities in a 
Charge- Transfer Insulator 


It has already been discussed that electron spectroscopy experiments suggest that the 
parent compounds of the high-T c superconductors are CTI’s rather than MH insulators. This 
itself makes the study of static impurities in a CTI interesting. Otherwise, in view of the 
experiments on static, nonmagnetic-impurity-doped systems, although some of their proper- 
ties have been understood within a one-band model, as discussed in the previous chapters, 
others require a more detailed picture like the three-band model for their understanding. 

In the one-band Hubbard model in hole picture, the defect state is formed slightly below 
the upper Hubbard band (UHB), in the Hubbard gap. This being the lowest unoccupied 
state for the holes, the first added hole goes into this state and the moment is lost. Due to 
thermal excitation of the hole from the defect state to the UHB at temperatures comparable 
to the energy difference between the two, and also due to quantum fluctuations, part of the 
moment can be regained as discussed in Chapter 5. But how, for any kind of nonmagnetic 
impurity, a moment of ~ 1.2 \jlq per impurity is obtained irrespective of the doping level even 
at low temperature, and why this moment is ~ 35% less than the theoretically expected 
value [58] - these questions cannot be answered within a one-band model. 



Static Nonmagnetic Impurities in a Charge-Transfer Insulator 


In this chapter we consider the full three-band Hubbard model description of the Cu 3d 
;s in the copper-oxide planes in the insulating state and present a self-consistent numerical 
study of the system. In the first section we discuss our results for the undoped system. 
;he next section a study of the impurity-doped system is presented. 


1 Pure system 

We consider the three-band Hubbard Hamiltonian given by Equation 1.2 for 3d holes 
;he Cu0 2 planes. In this section we shall be talking about a half-filled system when we 
^ exactly one hole per copper site. 

For all our calculations we take tp d — 1.3 eV, U p = 0 eV, U pd = 0 eV, and e p = 3.0 eV 
1 e d = 0 eV so that the bare charge-transfer gap A = 3.0 eV, and do a numerical HF 
mlation on a a 10 x 10 Cu0 2 lattice using periodic boundary conditions in both x- and 
irections. A more detailed discussion on the choice of the values of parameters is given 
;he section on doped system. The numerical HF scheme is same as that for a one-band 
bbard model discussed in Chapter 5 the only difference being that now we have to take 
e of oxygen and copper sites separately. For AF symmetry breaking along the ^-direction 
y, the interaction terms can again be written as U + U (n isL )n^. In the site basis, 
trix elements of the Hamiltonian for spin a are given by, (i | H a | i) = e d + Ud{n ig ) if i is 
!u site and (i | H a | i) = e p + U v (n i5 ) if i is an oxygen site, (i \ H a | j) = ±t pd if i and j 
NN Cu and O sites and {l\ H a \m) — ±tpp if l and m are two NN O sites. After this, 
numerical calculation proceeds exactly as before till self-consistency is achieved. 

Exact numerical (Quantum Monte Carlo) calculations reveal that the three-band Hub- 
d model with a finite positive U d possesses an AF ground-state [71], In order to study the 
gnetic behavior of the system with changing U d we have calculated the sublattice mag- 
isation as a function of U d - Figure 6.1 shows the variation of sublattice magnetisation as 
unction of U d for different values of tpp. For t pp = 0 the system goes into an AF state 



6.1 Pure system 


63 



U d 


Figure 6.1: Sublattice magnetisation vs. U d for different t pp (all energies in eV). 

for any finite positive U d . This is because of nesting of the Fermi surface. When tpp ^ 0 
the system loses its nesting property and a critical value of Ud is required for the system to 
become antiferromagnetic. This critical value of U d increases with t pp . 

A study of the energy spectrum of the system shows that for large values of U d ( U d A) 
the system is a CTI. For an Lx L CuCA lattice, having L 2 Cu and 2 L 2 O sites, the first L 2 /2 
states form the LHB. The oxygen 2p band in the middle has 2 L 2 states and the UHB has 
another L 2 / 2 states. The 2p band is further split into three sub-bands similar to Zhang-Rice 
singlet-triplet splitting. The singlet sub-band lies below the bare oxygen energy, there are 
a group of non-bonding (NB) states at the bare oxygen energy e p and the triplet sub-band 
lies above these NB states. For U d < A but above the critical value so that the system 
is an AF, the UHB goes below the oxygen band. The charge-gap is maintained and the 
system is a MH insulator. Below the critical value of U d when the sublattice magnetisation 
vanishes, the charge-gap also disappears and the system becomes a d-type metal. There are 
no singlet-triplet splitting in the oxygen band in these cases (Figure 6.2). The situation in 



a) 


NB 



jure 6.2: Single-particle energy-bands for the three-band Hubbard model in the hole pic- 
e: (a) metal, (b) Mott-Hubbard insulator, (c) charge-transfer insulator with a singlet- 
jlet splitting. Shaded regions denote the occupied states. [N](A)B=[non](anti)bonding, 
singlet, T=triplet, E C t = renormalised charge-transfer gap. 
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intermediate range of values of U d is more complicated because the UHB and the oxygen 
band are strongly mixed. 

In the one-band Hubbard model and its strong-coupling analogue, the Heisenberg model, 
the nature of spin excitations has been studied extensively. It is well known that in these 
systems the spin-wave spectrum is gapless with the existence of a Goldstone mode and the 
dispersion for the long-wavelength spin-wave modes is linear in the momentum Q. As already 
discussed, the spin-wave modes are obtained from the poles of the time-ordered transverse 
spin propagator (T[S~(i,t),S + (j,t')]) which, in the random phase approximation, has the 
form [x -+ (fi)] = i_yj|^b(n)] • The matrix elements of [x 0 (fl)j, the bare antiparallel-spin 
particle-hole propagator, are given by, 


/ OO rjr j 

- «). (6-i) 

-oo Zj i 

where g ’ s are the Green’s functions in the AF state. In terms of the eigensolutions Ei, (pi in 
the self-consistent state, x°(0) has the form: 


. B -£ Er 4,4 

lx ( )Jii ~ E„-E ml -(i + £ Er nu-Ev + n- ( 

Since only the Cu spins form AF order, spin-waves will have amplitude only on the Cu sites, 
and so we calculate [x°(^)]r? connecting only the Cu sites. In Equation 6.2 the site indices 
i and j thus run only through the Cu sites, but the eigenvalue indices l and m run through 
all three bands. 

All information regarding the nature of the spin-wave excitations is again contained in 
the eigensolutions of [x°(^)]> 35 the full RPA susceptibility can be expanded in terms of the 
eigenvalues A(Jf2) and the eigenvectors | <p\(Q.)) of [x°(^)] as: 


Dr + (ft)] = £ 

A 


i 

1 - E d A(Q) 


(6.3) 
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3 spin- wave energies are obtained from the pole 1 — ?7dA(f2) = 0 in the RPA spin suscep- 
iity. and the eigenvector yields the spin-wave amplitude. 

In the three-band model, as in the one-band model, the largest eigenvalue A max (G = 0) 
the [x°m matrix for Q = 0 is found to be exactly equal to confirming that the 
n-wave excitations are gapless. Nature of spin-wave modes can be studied by obtaining 
; angles of rotation ft of local spin vectors from the spin-wave amplitudes fa and the local 
gnetisation SI using [3]: 

ft = sin -1 (6.4) 

e flo — 0 mode corresponds to a spin-wave excitation which has equal amplitudes on 
sites and rotates all the spins by equal angles and is identified as the Goldstone mode, 
us, as in the one-band model, the spin-wave excitation in a three-band model is gapless 

1 there exists a Goldstone mode as one one would expect in a system with continuous 
n-rotational symmetry. 

2 Impurity- doped system 

We have studied the problem of static nonmagnetic impurities doped in the three-band 
bbard model. As in the one-band case, we model the impurities by a high on-site repulsive 
.ential V to restrict the Cu 3d holes from occupying that site. The Hamiltonian for the 
)urity-doped system is, 


H impure + V a \cr a I(n (6.5) 

/,cr 

216 impure is the Hamiltonian for the pure system given by Equation 1 . 2 , and I runs over 
impurity sites. 

The values of most of the parameters of the model are related to electron spectroscopy 
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and are obtained from comparison with the results of these experiments. Different kinds of 
theoretical techniques give slightly different sets of values. But a generally acceptable set of 
values is [46] t pd = 1.3 eV, t pp = 0.65 eV, U d = 8.8 eV, U p = 0 eV, U pd = 0 eV, A = 3.0 
eV. In this section we work with this set of values of the parameters. The pure system at 
half-filling is a CTI with this set of par am eters. 

Starting from the pure system at half-filling one Cu 3d hole is removed for each added 
impurity. For each impurity, depending on the sublattice of its substitution, a self-consistent 
solution of this Hamiltonian shows that in the spectrum of one particular spin there are 
three states outside the bands and for the other there is one outside. Consider an impurity 
on an A-sublattice site (having majority of f-spin holes) in an L x L system. In the f-spin 
spectrum there are (L 2 / 2 — 1) states in the LHB. One state goes out to the impurity state 
which is the highest in energy (~ V). The oxygen band and the UHB have (2 L 2 — 1) and 
(L 2 / 2 — 1) states respectively, one state from each coming down to form two defect states. 
The L 2 /2-th and the 2L 2 -th states are the two defect states. In the spectrum for the F-spin 
hole, the LHB and the oxygen band remain unaffected. The UHB has (L 2 / 2 — 1) states, 
one state going to the impurity state. The situation is simply reversed if the impurity is 
placed on a B-sublattice site. The energy-spectra for the two spins for a single impurity on 
an A-sublattice site is shown in Figure 6.3. 


Coming to the nature of these states, we find that the impurity state is essentially 
site-localised on the impurity site, whereas the lower defect state, coming out of the oxygen 
band, has amplitude mainly on the oxygen sites around the impurity, and the upper defect 
state has amplitude mainly on the four Cu sites neighbouring the impurity. In this respect 
the upper defect state is similar to the defect state in a one-band model. This state, when 
occupied by electrons, creates a local moment residing predominantly on the four Cu sites 
neighbouring the impurity. With only a nonmagnetic impurity on an A-sublattice and no 
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extra added hole, the system has (L 2 / 2 - 1) f-spin and L 2 j 2 4-spin holes. So the LHB is 
completely filled and all the other states are empty of holes. The lower defect state being 
the lowest unoccupied state for the holes, the first added f-spin hole goes into this state and 
subsequent holes go into the oxygen band. So even at finite concentration of hole doping 
the upper defect state remains occupied by electrons. Thus unlike in one-band model, the 
local moment in three-band model is robust with respect to hole doping. In Figure 6.4 we 
show the variations of the two energy gaps — between lower defect-state and oxygen band 
and the upper defect-state and the UHB. It is interesting to note that although the lower 
gap varies considerably, the upper gap is relatively insensitive to the variation of V. 

The magnitude of the local moment is proportional to the density of the defect state on 
the copper sites. We find that the total density of the upper defect state on the copper sites 
is less than 1, and that for a larger value of t pd the reduction is larger. In our calculations, 
for t pd = 1.3 eV about 21% of the defect state density goes out to the oxygen sites, whereas 
for t pd = 0.75 eV, only 7% of the density goes out. This supports the view that part of the 
local moment is transferred to the oxygen sites because of Cu-0 hybridisation, which leads 
to the observed moment being less than the theoretically expected value. Figure 6.5 and 
Figure 6.6 show total density of the upper defect state on Cu sites for t pd = 1.3 eV and 0.75 
eV respectively. 

The density is larger for smaller t pd but is interestingly found to be essentially constant 
over the range of V studied in our calculations. This leads to the significant conclusion that 
the defect-state-induced moment is more or less independent of the details of the nonmagnetic 
impurity within our model of the impurity potential term and leads to an understanding 
of the observation that independent of their detailed electronic nature, all nonmagnetic 
impurities, Zn +2 , Al +3 , Ga +3 , give moments of same magnitude. 

We have also studied the nature of spin-wave excitations in the impurity-doped system. 
We find that the gapless Goldstone mode associated with spontaneous symmetry-breaking in 
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the AF ground state is preserved in presence of impurities and even when the concentration 
of impurities on the two sublattices are not equal. This is because the continuous spin- 
rotational symmetry of the Hamiltonian, which is broken in the AF ground state, is still 
preserved. Wan et. al. [5] have obtained two modes of spin- wave excitation in case of 
unequal concentration of impurities on the two sublattices. We believe that the gapless 
mode in our calculation corresponds to their acoustic mode. 



hapter 7 


tatic Magnetic Impurities in a 
lott-Hubbard Antiferromagnet 


The effects of nonmagnetic impurities on a Mott-Hubbard AF was studied in Chapters 2 
1 4. The problem of magnetic impurities in an AF has been studied within a Heisenberg 
del by Wan et al. [5]. But to our knowledge, no such study of magnetic impurities 
i Mott-Hubbaxd AF exists. In this chapter we develop several formalisms, representing 
erent situations, to treat magnetic impurities in a Mott-Hubbard AF. 

In the first section we look at the problem of a single-orbital magnetic impurity in a 
gle-orbital Hubbard model. The impurity is represented by a modified hopping between 
impurity-orbital and its nearest neighbours, and its effects on the magnon-spectrum is 
died both within a perturbative and an exact-eigenstates analysis. In the second section 
extend the spin-independent potential used in the case of a nonmagnetic impurity to 
pin-dependent potential for a magnetic impurity and study its effects on the electronic 
ictra and spin-wave properties. In the last section we study the effects of higher-spin 
mrities on a spin 1/2 AF by multiple-orbital representation of the impurity site. We 
dy their effects on electronic properties and on magnon properties through a numerical 
' study on a finite-sized system. 
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7.1 Single-orbital magnetic impurity 

In this section we consider a single-orbital magnetic impurity in an AF host which is 
described by the Hubbard model on a square lattice with one orbital per site at half filling. 
The host Hamiltonian is again given by: 

Ho t ^ T Q’jcr&iar') "b U ^ * Tltf (7-1) 

{ij)<r i 

where t is the NN hopping strength and U the on-site Coulomb repulsion. 

We model the single-orbital impurity by taking the hopping term if between the impurity 
orbital and its NN orbitals to be different from the NN hopping terms t in the host system. 
The Hamiltonian with such an impurity can be written as below, where the sum is over all 
the impurity sites I and their nearest neighbours J, and 5t = t! — t. 


H = Ho + St J2 (“U + «U) (7-2) 

(JJ)cr 

In the strong-coupling limit the Mott-Hubbard AF with //-orbitals per site maps onto 
an S = A/"/2 quantum Heisenberg AF with NN coupling J = qy-. So the present situation 
represents an S' = 1/2 impurity in an S = 1/2 host AF with modified couplings J' = ~ 
between impurity-spin and its NN’s. We treat the host AF in the HF approximation as 
discussed in Section 2.1, and treat the effects of magnetic impurities on this HF state. 

7.1.1 Perturbative Expansion 

We first consider the perturbative technique which is exactly similar to the one used in 
the case of a nonmagnetic impurity discussed in Chapter 4. The idea is to calculate the cor- 
rection [t5y^] to the [x^j matrix in powers of 8tjt , from that to calculate the correction to the 
eigenvalues of [x°] perturbatively and hence to calculate the renormalised magnon energies. 
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Diagrams which contribute to [5xu] and t0 d rst or d er in & are shown in Figure 7.1, 

where the continuous lines represent the host Green’s functions and the dashed lines are the 
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Figure 7.1: 

hopping perturbations 6t between the impurity site I and the nearest neighbour J. The 
upper and lower lines are respectively the f-spin and 4--spin propagators with frequencies oj 
and u — ft. 

We now discuss the evaluation of the matrix-elements of the [hx°(ft)j matrix in the 
strong-coupling limit wherein, as before, terms only upto order t 2 / A 3 are retained. The first 
term in <5x?j is given by the expression, 

if — \ - .5 1. - 1- 

J 2?r \u) — A - — iri uj — A + + irjj \u) — ft — A~ — ir) a; — ft — A+ + irj 

where we have taken terms upto 0(1) in the local Green’s function and upto 0(t/ A) in 

the NN Green’s function, and Ef = A* = ±A. The integral can be easily evaluated 

using contour-integral method. Adding equal contributions from the two diagrams we have, 
rr..(n _ l St t 2 

l d X J/J - 

To calculate {6x°n\ we consider 8G , the perturbation-induced correction to the Green’s 
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function. As discussed in Section 4.1 


(0)]y — i J (dui/2'ir) — fi) + g^^SG^iu — f2)j (7-4) 

Since x° is the antiparallel-spin, particle-hole propagator, we need to examine the change in 
spectral weights of Green’s functions in the lower and upper Hubbard bands. Therefore we 
consider first the correction to the local spin-up Green’s function to first order in hopping 
perturbation, 5G* : = 2J2j gJjStgjj, where the sum is over all 2 nearest neighbours. We take 
the impurity site I to be on the A sublattice. In the strong-coupling limit, again substituting 
the appropriate expressions upto order t / A for host Green’s function matrix elements, we 
obtain: 





—t/2A tj_ 2A 

U — A - — ig a ; — A + + if] 


(7.5) 


Now the change in spectral weight in the lower Hubbard band is obtained by integrating the 
imaginary part of the Green’s function over the lower Hubbard band. Contribution comes 
only from the particle-hole possibility, the second-order pole at A~ yielding no contribution, 
and we obtain: 





(7.6) 


From this change in the spectral weight it is straightforward to obtain the correction 
5x°n to the particle-hole propagator. Since the 4.-spin spectral-weight in the upper-band is 
1 to leading order, these give the correction [5x°]// = "~f -A- similar calculation gives 

[Sx°]jj = ~ 2 Ta J* So ’ u P to 0(j?)> the reievant matrix elements of the [$x°] matrix are: 


M// = - 


z 5t t 2 

27A3’ 


[5x%j = Mji = jj = - 


1 5t t 2 

27A3‘ 


(7.7) 
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identical result is obtained if the impurity is placed on the B-sublattice. 

As we have already discussed in Section 4.1, in the long-wavelength limit, the eigenso- 
.utions of the [x°] matrix are: 



The first order correction to the eigenvalue is given by 


«' = {<? | fe°] I Q) = i (a|<Sx?, + a^zi^Xu + Aja^T^Xj, + fyzSx] j) ■ (7.10) 

If the concentration of impurities on the A- and B-sublattices are x A and xb respectively, 
then summing over the contribution from all impurities, we have, 


(S\ l ) = ~(x A +x B )z 


5t t 2 

7 A 1 


1 (—) 91 

2 V2J/ + 4 


(7.11) 


The renormalised magnon energy is now obtained from the solution of the equation 1 — 
U{ A 0 + SX 1 ) = 0, and upto first order in the total impurity concentration 2x = x A + x B , we 
obtain: 


1 + 2 z.x . — 
t 

where q = V2JQ is the magnon energy for the two-dimensional host antiferromagnet in 
the long-wavelength limit. The result agrees exactly with the calculations on the Heisenberg 
model [5] in that there are no singular corrections to the magnon energy in the case S' = S , 
and the correction is proportional to St/t = (1/2)5 J/J. In Reference [5] the magnon energy in 
the long-wavelength limit was obtained as: e Q (x) = eg(0)[l -I- 2 x(j - 1 )/j)], where j = J'/J. 



CIq = 
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7.1.2 Exact-eigenst at es analysis 

In the exact eigenstates method, we do a numerical HF study on a finite-sized Hubbard 
model as discussed in Chapter 5. For the case S' = S = \ which we consider here, the 
magnetic impurities are represented by a modified local hopping term, i.e., (i | H | j) = t 
if none of i and j is the impurity site, but = f if one of them is the impurity site and the 
other is a NN. Once the self-consistent AF state is obtained, the eigensolutions are used to 
construct the [x°(C)] matrix as discussed in Chapter 6 and the magnon modes are obtained 
from the the eigensolutions of the [x°] matrix. The energy Q, n on the n-th magnon mode 
is obtained by solving, 1 — U\ n (fl n ) = 0 for the appropriate (n-th from the top) eigenvalue 
A n (f2). The root of the equation, A n (f2 n ) = 1/U is determined by obtaining A n (f2) for closely 
spaced values of on both sides of the root, and finally linearly interpolating between them. 
Suppose (A*,f2*) and (A^,f2^) are two sets of values for two energies very close to, and on 
either side of the root, then the root Q n is determined from, 

A(n„) = i = Ai + A^A(n„-n 1 ) (7.i3) 

Here we focus on the short wavelength (high O) modes. The magnon spectrum is shown 
in Figure 7.2 for different values of St/t for a 10 x 10 lattice with U/t — 10, and clearly 
shows that precisely one magnon state at the upper end of the spectrum is split off from 
the magnon-energy band. This is the short-wavelength mode corresponding to local spin 
deviation at the impurity site. Furthermore, it is seen that the energy difference of the 
split-off state from the upper end of the spectrum increases roughly in proportion to St for 
small St’s. This is easily understood as making a local spin deviation on the impurity site 
costs an energy Ae oc {J 1 — J ), i.e., Ae oc St/t to first order in St/t. The localised nature of 
the split-off mode is shown in Figure 7.3 where the square of the magnon wave function for 
this mode is plotted on the lattice grid. 
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Figure 7.2: Magnon spectrum for a single-orbital magnetic impurity. 
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Figure 7.3: Magnon wave function (lower) and sublattice magnetisation (upper surface). 

7.2 Spin-dependent impurity potential 

Nonmagnetic impurities such as Zn, Al, Ga, etc. in the Mott-Hubbard AF were rep- 
resented via spin-independent impurity potential term V Y,ia a \<r a io hi the Hamiltonian. A 
natural extension for magnetic impurities is in terms of a spin-dependent impurity potential 
term, and we therefore consider the following impurity term in the Hamiltonian 

H£S = P- 14 > 

I 

where T/ = (a/f an). A spin-independent impurity potential eo can be included for gener- 
ality, however, we shall consider the limit V > e 0 , so that the potential for spin a fermion is 
V a ss —oV. We choose V to be positive if the impurity is on an A-sublattice site, so that 
is very low and V t is very high. The sign of V is reversed for impurities on B-sublattice sites. 
This choice of potential ensures that the magnetisation on the impurity sites follows the host 
AF ordering. Such a spin-dependent impurity potential can arise from a coupling — a.S- imp 
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between the itinerant fermion spin a and the static magnetic impurity spin -Simp. Within a 
classical picture the impurity spin provides a local magnetic field at the impurity site. Since 
experiments on high-T c cuprates doped with static impurities show the impurity spin to 
be antiferromagnetically coupled with the host spins [58], we take the local field direction 
to be along the local magnetisation direction (z), and this yields the -a z V potential. We 
thus have an effective one-orbital-per-site model for the magnetic impurity in the AF host in 
which the spin-f (-4-) fermions feel a strong attractive (repulsive) potential at the impurity 
site. 

The remaining analysis of this section closely follows the earlier study with nonmagnetic 
impurities [64]. Within a T-matrix approach the Green’s function for the two spins can be 
written as 


Energies of impurity-induced states are then obtained from the singularities in the T-matrix, 
and are given by solutions of gjj(ui) = 1/V4- The local host Green’s function in the AF state 
was obtained as gjj(uj) = jj — o' A) / (w 2 — E~) for I 6 A where a = ±1 for f and 4-spins 
respectively [63]. 

A plot of gJj(cj) vs. uj is shown in Figures 7.4 and 7.5 for the two spin cases together 
with 1/V a . Intersection of gJ/(uj) with the —1/V line gives a impurity state for spin-t at 
an energy ~ -V. Similarly, there is an impurity state for spin-4 from the intersection of 
5//( w ) and 1/V at an energy ~ V. For large \V\/U, these impurity states are essentially 
site localised on the impurity site, and therefore decoupled from the system. Similarly when 
the magnetic impurity is on a B-sublattice site, again no defect states are formed in the 
gap, only impurity states are formed. This follows immediately from the spin-sublattice 
symmetry which results in gfj(uj)\i^B — gn( u )\i&A- Also the signs of the impurity potentials 
for the two spins are reversed, as mentioned earlier. As a consequence the intersections with 
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l[V a yielding the two impurity states are at the same energies. 

Thus a significant difference for the magnetic impurity case, when the impurity spin 
is antiferromagnetically coupled to the neighbouring host spins, is that there are no defect 
states formed in the Hubbard gap, rather there are only impurity states formed far removed 
in energy from the Hubbard bands. 

For magnetic-impurity doping in the AF host the fermion number is unchanged, unlike 
the case of nonmagnetic impurities where one fermion is removed for every added impurity. 
However, in the limit V — > oo the impurity site gets decoupled from the host in much the 
same way, so that an identical impurity-induced perturbation is obtained. This 

is because while for the nonmagnetic impurity placed on, say, an A-sublattice site, it was 
the removal of one spin-t fermion which reduced the minority spin-t density by t~ / 4 A 2 on 
the NN B-sublattice sites, for the magnetic impurity case, it is the strong attraction due to 
potential — V at the impurity site which leads to this depletion. 

We now turn to the renormalisation of magnon energy. In the limit of V -4 oo, 
all extended-state wave functions have vanishing amplitude on the impurity site, and the 
impurity-localised states at infinite energy have vanishing contribution to the particle-hole 
process for small (order J) energies. Therefore, [x°(f2)]y must vanish if i/j = I, so that, 

[tx°(n)h = -[Xhost(^)3o (7.16) 

The only other terms in [6x°] to O(^) are 6x°jj where J’s are NN to I. The small 4-spin 
particle density of the order of ^ sitting on the A-sublattice site I is repelled due to high 
potential V to the four NN sites J, so that the change in each one of them is 
Similarly the small f-spin hole density of is also removed to the neighbouring sites so that 
^ n ft =: 4 a 2 - corresponding to the shift in 5nj t to the impurity site. The situation is exactly 
similar to the case of nonmagnetic impurity except that we are now working in the electron 
picture. Proceeding as earlier (4), we have 5xjj = , and so the various matrix elements 
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of the correction to the bare, antiparallel-spin particle-hole propagator are as follows: 


[ s x‘<- l i\ n - -JJ + 2^11 7- jj) (7.17) 

(7.i8) 

These are identical to the case of a nonmagnetic impurity in an AF and lead to an 
identical result for the renormalised magnon energy. For low impurity concentrations and 
in the long-wavelength limit, we again have a logarithmically divergent reduction to the 
spin- wave energy in 2D, 


= 


V2JQ 
= V2JQ 


1 — X — 4x J2 — 2 — 7 ~r 
N ^ Q' 2 - Q2 

1 - X - Ax (— In ^ ^ 

U Q 27 


(7.19) 


where Q c is an upper momentum cutoff of order 1. In two dimensions there is a logarithmic 
singularity as Q -4 0, indicating a logarithmically divergent reduction to the magnon velocity. 
Impurity scattering also introduces a damping of spin-wave modes. 

Therefore we conclude that irrespective of whether the impurity is magnetic or nonmag- 
netic, if a site is decoupled from the system by putting a large on-site potential - attractive 
or repulsive - it acts as a ‘strong’ perturbation to the spin-wave excitations and leads to sin- 
gular corrections to spin-wave energy in 2D. The hopping-perturbation, on the other hand, is 
a ‘weak’ perturbation for long-wavelength modes and only causes a momentum-independent 
shift in the spin-wave energy proportional to the impurity concentration. 
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7.3 Multiple-orbital representation of a magnetic im- 
purity 

We can represent higher-spin magnetic impurities (S' > 1/2) within the Hubbard model 
by incorporating multiple orbitals at the impurity site. An S' = y impurity can be repre- 
sented by putting N orbitals on the that site. A HF numerical study can be done on this 
system with the only difference that we now have an enlarged basis including the N impurity 
orbitals. We c an reach a self-consistent state in this enlarged basis and from that calculate 
the magnon energies exactly as in Section 7.1. We focus on Goldstone mode, electronic 
feature of the magnetic impurity and magnon-energy renormalisation due to scattering off 
the impurities. 

As the simplest approximation, we assume all the impurity orbitals to be connected to 
the NN’s through the same hopping f, and to involve the same on-site Hubbard repulsion U 
as in the host system. The Hamiltonian for such a system can be written as: 

/ / 

H t (cL^Ujer h.c.) + U y ] t ^ ) (^laa^-Jv "t" h.C.) + U y ) (7.20) 

{ ij)cr i (IJ)aa la 

Here J are the NN sites of the impurity site J, the primes represent sums over sites except 
the impurity sites, and a = 1, 2, ... N denote the N orbitals at the impurity sites. Since the 
N impurity orbitals are independent, in the strong-coupling limit this model can be mapped 
onto a quantum Heisenberg model with the Hamiltonian: 

Heff = Jj: Si-Si + j £ sUs} + SI + - - - sf) (7.21) 

(ij) (H) 

where, as before, J = and Sf = t /ft (a/ 2)t /a with a = 1,2,.. ., JV axe the impurity 
spin operators constructed from the fermion operators. The “impurity spin” S^ p = S} + 
Sf -j f- Sf is not exclusively a spin N/2 object. However, in the HF AF state, since all 
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the impurity orbitals are connected to the same set of NN’s through the same hopping, they 
will preferably be occupied with electrons of the same spin, thus effectively projecting out 
possibilities of other spin configurations. 

We expect the Goldstone mode in the transverse spin fluctuations to be retained even in 
presence of impurities as the continuous spin-rotational symmetry is preserved as is evident 
from the effective Heisenberg Hamiltonian in Equation 7.21. One way to understand trans- 
verse spin fluctuations in the Hubbard model is in terms of the rotation angles = sin -1 ff, 
where and s\ are the magnon-mode amplitude and the local magnetisation at the site i 
[70]. Now, in order to treat the magnetic impurity as a single object, we simply add the mag- 
netisations and the magnon amplitudes on all the impurity orbitals to obtain sj = Ea=i s ia 
and 0/ = J2a~i ®r«- In our numerical HF studies we find that for 0 = 0, the largest eigen- 
value of the [x°(0)]-matrix is exactly equal to jj and the rotation angles for this mode are 
identical on all sites, including the impurity sites with 0/ = sin -1 for N = 2, and 3, thus 
explicitly confirming the existence of gapless Goldstone mode. 

7.3.1 Electronic features of multiple-orbital magnetic impurity 

We now consider some interesting features in the electronic spectra in presence of the 
multiple-orbital magnetic impurities we are considering. For concreteness, we consider the 
impurity to be on an A-sublattice site. In the HF state in the strong-coupling limit, the 
f-spin and .[.-spin amplitudes are 1 — -Jy and respectively on an A-sublattice site. This 
continues to hold even for the impurity orbitals as they are connected to four NN’s through 
hopping f , and a density of is transferred from each one of them to each one of the NN s. 
However, in our numerical HF calculations we find that when y is not too large there is an 
enhancement of local magnetisation on the impurity orbitals. For example, for three orbitals 
on the impurity site, (i.e., S' = 3/2) the magnetisation is uniform for a- = 25, but there is 
an enhanced magnetisation on the impurity sites when ~ = 10. This is because the hopping 
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of an f-spin electron from one of the impurity orbitals to the NN s suppresses the hoppings 
of the electrons in the other impurity orbitals and this is a higher-order (in powers of jj) 
effect. So the "["-spin electrons in the impurity orbitals tend to get localised and hence get 
partially decoupled from the system. 

Within HF approximation this enhanced magnetisation on the impurity orbitals leads 
to a lowered HF potential for f-spin electrons on these orbitals. This relatively lower HF 
potential can be treated as a weak attractive potential K// = U(ri\ — n£ ost ) for the f-spins 
on the impurity site and it is easy to see from our earlier T-matrix analysis (??) that this 
will lead to a localised state just below the LHB by solving the equation 1 — V e ffg]j(u>) = 0. 

There is another mechanism in case of multiple-orbital magnetic impurities which leads 
to localised states on the impurity site. Since all the impurity orbitals are connected to the 
four NN’s through the same hopping t, any localised state on the impurity site such that its 
amplitudes on the impurity orbitals add up to zero, i.e., the hopping contributions from the 
different impurity orbitals cancel exactly, will be an exact eigenstate of the HF Hamiltonian. 
This is energetically made possible by lowering of the HF potential for f-spins at each of 
the impurity orbitals. So, for an A r -orbital magnetic impurity, a site localised state on the 
impurity site with amplitudes fa, fa, • ■ ■ , 4>n on the N impurity orbitals will be an exact 
eigenstate of the HF Hamiltonian if the following condition is satisfied: 


4 >\ + fa H ¥((>n = 0. (7.22) 

Actually there are two conditions involving these N quantities. One is given by Equa- 
tion 7.22, and the other is the demand that the site-localised state be normalised, i.e., 

|dh| 2 + \fa\ 2 4 \4>n\ 2 — I- (7.23) 

For a spin 1 impurity with N = 2, there is a unique solution with fa = -fa = 
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satisfying these two conditions. So, for a two-orbital magnetic impurity, a site-localised state 
which is antisymmetric between the impurity orbitals will be an exact eigenstate. However, 
for N > 2 there are infinite possibilities for the <p Q ' s. Any suitable linear combination of these 
infinite possibilities will be a site-localised exact eigenstate leading to partial decoupling of 
the impurity site. 

7.3.2 Magnon-energy renormalisation 

Now we focus on magnon-energy renormalisation due to scattering off the impurities. 
We have seen in the case of a nonmagnetic impurity and a magnetic impurity represented by 
a spin-dependent potential that decoupling of a site from the system causes strong scattering 
of the spin waves and leads to a diverging correction to the spin- wave velocity in 2D. It will 
be interesting to see how the spin-wave velocity gets renormalised in this case. In their work 
on impurities in a Heisenberg AF, Wan et al. [5] found a logarithmically diverging correction 
to the spin-wave velocity in the long-wavelength limit in 2D for S' S. So it is of interest 
to see the scaling of spin- wave velocity renormalisation due to scattering off the impurities 
with system size. 

Spin- wave modes in a finite-sized system can be assigned wave-numbers Q = x + n l 
[70]. The first mode after the Goldstone mode corresponds to n x = 1, n y = 0 or equivalently 
n x = 0 and n y — 1 so that Q = Similarly, the second mode corresponds to Q = qy\/2- In 
Figure 7.6 we plot the spin- wave energy for these modes for a 12x12 system with two spin 1 
impurities, represented by putting two orbitals on each of the impurity sites. The spin-wave 
energy is found to be linear in Q for small Q’s even in presence of impurities. 

As the spin- wave energy is proportional to the wave-vector, the spin- wave velocity v 
can be defined by = vQ in the long-wavelength limit. The spin-wave velocity has been 
expressed in the form v(x ) = w(0)(l — ax) [5], in the low-impurity-concentration limit, a 
has two parts and can be written as a = ao + an, where a 0 is a constant independent of 
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Figure 7.6: Spin-wave energy for small Q’s in presence of impurities. 



x ( in % ) 

Figure 7.7: Renormalisation to the spin-wave velocity as a function of impurity concentration 
for different system sizes — 8x8 (diamond), 12x12 (plus) and 16xl6(square). The points 
are the calculated values and the lines represent the best fits. 
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Figure 7.8: Plot of a vs. lnL. The points are the calculated values and the line represents 
the best fit. 

the system-size and oll depends on the system-size L. If there is a logarithmic divergence in 
the renormalisation to the spin-wave velocity in the long-wavelength limit, then we should 
see an increase in the absolute value of the slope of the spin-wave-velocity-renormalisation 
vs. impurity concentration curve and should have a cc lnL. In Figure 7.7 we show the 
renormalisation of spin-wave velocity with doping concentration for three different system 
sizes — 8x8, 12x12 and 16x16. It is clear that for the small impurity concentrations we have 
studied, the spin-wave velocity decreases linearly with increasing impurity concentration and 
the ^-dependence (absolute value of the slope) increases with system size. This implies a 
diverging correction in the limit of large system. In Figure 7.8 we plot a, the slope of the 
impurity concentration vs. spin- wave velocity curve, as a function of In L and find that the 
variation is linear, which signifies a log-divergence. 
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Appendix A 


Local Green’s function and position 
of the defect state 

For I e A, 

9ii(v) = 


Q t 2 W | Q t 2 (£) 

u> + Ej; yj — Ej: 
u ( af{k ) + af(k)) - % ( af(k ) - af(k )) 

" Y w2 -^I 

We have af 2 (k) + af 2 (k) = 1, and from Equation 2.3 we have af 2 (k) - af 2 (k) = Jk 
Using these we can write, 

*}»(") = Ejfla- (A-D 

£ £ 

Similarly, for the down-spin Green’s function, using a® (k) + a® (k) = 1, and a® ( k ) — 
af 2 (k) = A ; we can write, 



a® 2 (fc) | af(k) 
u-Ef u-Ef 

k k 
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Local Green’s function and position of the defect state 


9iiH = E 

k 


u) + A 


These two can be combined to write, 


fn = E 

k 

where a = ±1 for f- and 4-spins respectively. 
For the f-spin Green’s function, 


u — a A 

a- 2 -4’ 

k 


(A.2) 


(A.3) 


S//M — .. E 


u) — A 
N w 2 — E\ 

k k 


r* ef + (A 2 — A 2 ) 

k k ' 




lne de 


e 2 + (A 2 — a; 2 ) 


(A.4) 


where we assume a logarithmic density of states in 2D. Evaluating the integral we have, 


(a.5) 

As u — f A from below, this — 0 as \/A — w in (A — u). 

As u; -4 — A from above, this blows up as (A -f A) -1 / 2 In (A + A). 

Since a solution of the equation g]j(oj) = 1/V gives the position of the defect state, and 
the defect state is formed close to the upper Hubbard band in the limit of large V, from the 
first of the above relations we have, 


\/v. 


(A.6) 


where tj is the difference between the upper band-edge (A) and the defect state energy (Ed)- 
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Defect state wave function 


The defect state wave function on a site i is, by definition, = d] I (E d ) /yj — — j Z=Ej- 
The denominator will give just a constant. So, up to a multiplicative constant, the defect 
state wave function is simply gJj(E d ). 

As we have assumed, / G A. First we consider i GA. gh{Ed) is then given by the Fourier 
transform of g\ A [k,E d ). 


9ii(Ed) = E 


af(k) 


E d 


Ef 

k 


+ 


a® 2 (£) 


E d 


Ef 

k 


Ak.f 


(B.l) 


where f = fi - fj. 


In the strong-coupling limit af 2 (k) = 1 


et 

4A 2 


and a® 2 (£) 


4^2 , so that 


al, (Bt) = £ 

k 

We have E d = A-rj. With this, ~ 2A ( X “ &) and e^e f ~ 7§k ( X “ Tffiz) • 

if it » V- 

From these we have, 


tt 

4 A 2 


+ 


4A 2 


E d ~ Ef E d - Ef 


Ak.f 


(B.2) 
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Defect state wave function 


<?!/(£*) = ^E 


V 


2A 2 e 2 J 2 A 


Ak.f 


(B.3) 


As 17 -» Jj, in the limit V -* 00, the defect state wave function on the A-sublattice 
sites are 0(^3) or smaller. 

For i G B, the defect state wave function is given by the Fourier transform of g^ BA (Ed) 
and in the strong-coupling limit, 


= e 

k 


E d - Ef E d - Ef 



(B.4) 


Proceeding as before, we obtain. 


*&<*> = s? 

k 

This shows that in the strong-coupling limit and with V — >• oo, the defect state has 
vanishing amplitude on the A-sublattice sites. Whereas, to leading order, the defect state 
wave function on the B-sublattice sites are given by <^(f) ~ Yk ^ l&lk ' T ■, U P to a multiplicative 
constant. 

On the other hand, when V is not too large so that 77 is not too small, we have E = 

k 

- — L — p-- This can be written as (l - > if V » ft- ■ In this limit, 

e^Ze¥ ~ ~jr (l ~ 2A^)- Proceeding as before, for i € B, we have, 


2 A 


L z k 


e k 

2A 


A k.f 


(B- 5 ) 


sl,(E d ) 


E 

k 






2A(2A - 77 ) 2 At? 


Ak.r 


(B.6) 


Fourier transform of €jr ~ (cos k x + cos k y ) is nonzero only when r is a vector connecting 
N T N sites. This shows that when V is not too large, defect state will be localised only on the 
four NN sites. 
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Perturbation due to hopping 


Here we present the corrections to the energy states of a two-site system due to hopping 
t. The notation we use is the same as in Reference [72]. We have four states - two at 


-A with eigenvectors 


1 


0 


" 0 ' 

0 

and 

0 

and two at +A with eigenvectors 

1 

0 


cos (9/2) 

0 

0 


1 

i — 1 * 

to 


0 


and 


0 

0 


sin(0/2) 

— cos(#/2) 

For the states at -A the V-matrix will be a 2 x 2 one with V {j = (i\[t}\j), where |i) 
and |j) are the eigenvectors of the Hamiltonian (Equation 3.2) with eigenvalues —A and [f] 
is the hopping matrix in the two-site, two-spin basis given by, 



0 0 -t 0 

0 0 0 -t 

-t 0 0 0 

0 -t 0 0 


(C.l) 


A straightforward algebra gives, 



Perturbation due to hopping 


r 0 —t cos (0/2) (0.2) 

^ ^ — icos(0/2) 0 

if 1 ! 

The [^matrix has eigenvalues -A ± t cos(0/2) with normalised eigenvectors 75 [ 1 J and 

^ 1 . This shows that the two states at -A are split because of the hopping to 

-A ± fcos(0/2) and the weights are equally distributed to the states at these energies. 

A similar calculation shows the states at +A to be split to + A±t cos(0/2) with weights 
getting equally distributed to the split-up states. 
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Hamiltonian for 

a 

five-site cluster 


H = 
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